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PREFACE 



The Department of Aeronautics, Naval Postgraduate School 
offers several courses at different levels on the dynamics of 
flight vehicles. These courses include studies of fixed wing 
stability, rotary wing stability, and the stability of 
missiles. In addition, there are general courses that cover 
the dynamics of bodies in six degrees of freedom, including 
the development of the equations of motion. 

There is currently no signle resource available that 
covers the material for the airplane and missile courses, as 
well as the introductory dynamics classes, and, as a result, 
the instructors must make use of large numbers of locally 
prepared handouts and/or assign several textbooks, only 
portions of which will be used. 

In order to provide a solution to this problem, the 
following material has been developed for use in the dynamics 
classes in the Department of Aeronautics. Although this 
material is, in places, of greater depth than what would be 
covered in the courses, it is felt that the value of this 
material is enhanced by its degree of thoroughness. 
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CHAPTER 1 



DYNAMICS OF MOTION 

1.1 INTRODUCTION 

This chapter is intended to be a formal review of 
mathematical and physical concepts that will be used to 
develop the theory of dynamic stability of flight vehicles. 

The starting point is the basic laws of Newtonian 
Mechanics, with a useful application to spring-mass-damper 
systems. We shall find later how the motion of any 
conventional flight vehicle is similar to those systems. 

We then proceed to rigid body equations of motion, 
described in non-inertial reference frames. 

1.2 NEWTON'S LAWS OF MOTION 

In any system of units, basic quantities and units are 
related through Newton's laws of motion. 

NEWTON'S 1ST LAW OR THE LAW OF COORDINATE SYSTEMS: 

If there is no force acting on a particle, its velocity 
remains constant. 

This law is often called the law of inertia and applies 
to inertial frames of reference. 

NEWTON'S 2ND LAW OR THE LAW OF DYNAMICS: 

The time rate of change of the linear momentum of a 
particle, is proportional to the net force acting on it. 
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The vector quantity 



( 1 - 1 ) 

is defined as the linear momentum of the particle, where m(t) 
is the inertial mass of the particle and v(t) is its velocity 
relative to the reference frame in question. 

The 2nd law is mathematically expressed as 




dp tU 
dl 



( 1 - 2 ) 



or 

F„«i ^ L v [i) ] 

If the mass m remains constant, we can write: 



(1-3) 



■ i m i 

or 



(1-4) 




= — on a GO 
3c 



(1-5) 



where a(t) is the acceleration caused by the net applied 
force 

The proportionality factor is expressed as V$ c and its 
value depends on the system of units being used. 

In the English Engineering system, one pound force 
(1-lbf) will cause to one pound mass (1-lbm) , an acceleration 
of 32.174 ft/sec^ Thus substituting into Equation 1-5 



- ~z~ Obi n) 32.174 C + Vsec 2 ) 

oc. 



( 1 - 6 ) 
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Hence 



9c “ 32.17-4 Ibrr? ft/lbf sec 2 (1-7) 

A few other units used for velocity, are given below with 
their corresponding conversion factors: 

1 mile per hour (lmph) = 1.4666 ft/sec 
1 knot = 1 nautical mile per hour = 1.6878 ft/sec 
NEWTON'S 3RD LAW OR THE LAW OF MUTUAL FORCES: 

If one body exerts a force on a second body, the latter 
exerts an equal but opposite force on the first. 

This law allows us to define the gravity force, or the 
weight of a body, that has a given mass (m) . It is the force 
that results due to the earth's attraction. 

Besides the three fundamental laws of mechanics, the 
Newton's Universal Gravitational Law states that two inertial 
masses ml and m2, separated by a distance R are subject to an 
attraction force according to the equation: 



f* = 9c G 



om, wt z 

~w~ 



( 1 - 8 ) 



where G, is the gravitational constant, which in the EE 
system is G = S.52*-lcf" Ibf --ftVfbw 2 



25 

If now m t = M e = 1.32 X 10 lbm (mass of the earth), 
then near earth's surface where R = R<? = 2.09 X 10 ? ft 
(radius of the earth), we can write: 



W - 9c 



G Mg 



'YU 
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(1-9) 



or 

( 1 - 10 ) 

where g = g c GM,<:/R| = 32 .174 ft/sec 2 , called acceleration of 
gravity. 

By Equation 1-10 we can see that 1-lbm, subject to the 
earth's gravitational field, experiences a force 1-lbf and 
hence it accelerates with 32.174 ft/sec 2 , near earth's 
surface . 

We can also define another unit of mass, the slug, as the 
mass that experiences an acceleration 1-ft/sec 2 , if subject 
to a force of 1-lbf. 

By this definition 

1 slug = 32.174 lbm (1-11) 

Consequently, we can derive the density (p) , as the mass 
per unit volume usually measured in slugs/ft. 

( P> = yokmew (1-12) 

People in the field use the (g) , as a short hand notation 
for force or acceleration. A 5-g force simply means a force 
5 times the body's weight. 

1.3 ANGULAR MOMENTUM 

A special significance to the 2nd law is given by the 
angular momentum approach derived simply by crossing the 
radius vector 7 , on to the 2nd law. 
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(1-13) 



r x 



c if 
di 



= r x 



P -riel 



We define 



G -nei = r x 



(1-14) 



the net torque or moment of the force about the origin and 

h = rxp (1-15) 

the angular momentum about the origin. 

Differentiating Equation 1-15 with respect to time f we 

get 

dk dr -* -* dp 

— = — xp + r X - 77 - 

dt di r di (1-16) 

The first term is zero, because 



dr -» 

dT*? 



v x (wu ) = -m (uxu) = o 



(1-17) 




Figure 1-1 

Moment of Force About the Origin 



So we finally get 

dk _ - (1-18) 

dt ' G " et 
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which tells us that the time rate of change of the angular 
momentum of a particle, equals to the net torque or moment 
acting on it. 

If no torque is acting on a particle, its angular 
momentum is conserved, likewise its linear momentum is 
conserved, if no force is acting on it. 

A simple example 

For a pendulum that swings in a vertical plane, having a 
mass (m) , the acting forces are, its weight W = mg and the 
tension of the rod T. The rod is assumed to have negligible 
mass. 




Figure 1-2 
A Simple Pendulum 

If 9 is the angle between the rod and the vertical axis, 
s, is the arc length corresponding to the angle 9 and 1, is 
the length of the rod, a direct application to the 2nd law 
(Equation 1-5), gives: 
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-mgsin© = on s (1-19) 

and equilibrium condition in the direction of the rod, gives: 

- 001 cj cos© = T (1-20) 

Equation 1-19, can be written 

s + gsm0 =0 (1-21) 

If 9 is small, we can assume 



Q - — 
y e 



5=19 



( 1 - 22 ) 



and 



9 = sin© (1-23) 

With this assumption. Equation 1-21, can be written 

9 + 4 - 9 =0 

11 (1-24) 

Equation 1-24, is a linear differential equation, with 
constant coefficients, which describes the motion of the 
pendulum. 

The same result can be derived, if we use the angular 
momentum approach. 

The angular momentum of the pendulum, about the pivot, is 
simply (£), cross the linear momentum (ms) . The time rate of 
change of the angular momentum, has to be equal to the net 
external torque. 



— (£*ws) = --fonqsin© 

Since l is constant (d8/dt = 0) , we can write 

fins = - £ vng s i n 9 



(1-25) 



(1-26) 
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or 



*s + W 0 sm 0 -= O 



(1-27) 



or 




(1-28) 



The minus sign in Equations 1-19 and 1-25 , indicates that 



and the angular momentum in Equation 1-25, increase. 

1.4 UNDAMPED HARMONIC OSCILLATOR 

In this section and in the following section, we will 
derive the equations and the solutions to the spring mass 
damper systems. 

The undamped oscillator is simply a spring mass system 
without friction, shown in Figure 1-3. 

Let k, denote the spring stiffness, m, the mass of the 
body, x(t), the displacement of the body from the equilibrium 
position. Notice that the weight of the body stretches the 
spring permanently, so that x e is the equilibrium position, 
after the body has hanged to the spring. 



if the angle (9 decreases, the acceleration s in Equation 1-19 






Figure 1-3 



Spring Mass System 
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(1-29) 



The restoring force according to Hook's law, is 

F s = - kx 

Then according to the 2nd law 

on x = - k x 

° r •• k 

x + — x =0 

As an example we will take 
m = 18 lbm 
k =40.5 lbf/ft 
If we let 



(1-30) 

(1-31) 

(1-32) 

(1-33) 



« a - 

™ (1-34) 

where is the so called natural frequency of oscillations. 
Equation 1-31 becomes 

x + X =0 (1-35) 

This is a linear differential equation with constant 
coefficients. Its characteristic equation is clearly 

j+ui ~ 0 (1-36) 

which has characteristic roots 

S\ = - j (1-37) 

The most general solution of Equation 1-35 can be 
expressed as 



xttMAj 



+ (A.) 



(1-38) 
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where A f and A_ are constants depending on the initial 
conditions. 

Expanding the exponential terms of Equation 1-38, we 
obtain 



x(i) = Oi COSUJn i *■[$ Sinujflt 

where 


(1-39) 


x — A + + A_ 


(1-40) 


& 

D 

a 

ii 

i 

>- 

i 


(1-41) 


If we now express 

smuj rt l = cos (iu„t + 90°) 


(1-42) 



we can give the following useful form to our solution 



M = COS ((jUy,-tf9 0 ) 


(1-43) 


where x wax is the maximum deviation from the 
position and 9« , the initial phase. 


equilibrium 



We can relate the quantities x^, a)t and , with the 
initial conditions, as follows: 

At time t=0, the initial displacement x„ , is given by 



Equation 1-43 . 

*o - cos ©o 


(1-44) 



and the initial velocity x 0 , is given by differentiating 
Equation 1-43 with respect to time and evaluating at t=0. 



*. = '’C/rtiax ^ n SI t) © 0 

As an example we will take 


(1-45) 


x 0 = 4 ft 


(1-46) 


x a = 32 ft/sec 


(1-47) 
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From Equations 1-44 and 1-45 , the initial phase and the 



maximum displacement can be found to be 



= ' U ( -(0,*. ) 


(1-48) 


, ( Z *o \*/2 




-wax * ~ ( X ° + ) 


(1-49) 



Note: In Equation 1-48, keep the sign, so that the 
quantity in brackets will give the right quadrant for the 
angle & . 

For our example 



' 3c s = 42 174 ’ - 79 39 2 

= W(^^-) = 456. 76’ = 2.59 r «A 

= ± 5.49 « 



(1-50) 

(1-51) 

(1-52) 



so that the response is 

x(i) - -5.49 cos ( 8.51 -t + 2.59 raol) (1—53) 

1.5 DAMPED HARMONIC OSCILLATOR 

Consider the spring mass damper system shown in 
Figure 1-4. Assume that the friction force exerted by the 
damper, is a linear function of the mass velocity. 
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xft) 1 
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k 



Uli 

o 
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r 
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=i c 



£ 



*1 



¥ 



Figure 1-4 

Spring Mass Damper System 



Let k the spring stiffness,- m the mass of the body, x(t) 
the displacement of the body from the equilibrium position 
x e , as it was defined in section (1.4) , and c, the 
coefficient of friction. 

The restoring force is again 

F s = - k x (1-54) 

and the friction force is 



- ~ C x 

Then according to the 2nd law 

on x =-cx-kx 



(1-55) 

(1-56) 



or 



c • 
X + — x 

-w 




(1-57) 
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As an example we will take 



If we let 



and 



18 lbm 


(1-58) 


40.5 lbf/ft 


(1-59) 


.2 . k 


*s« = 


(1-60) 


25 “» * s 


(1-61) 


so called damping coefficient, 


then the 


written 


x + 2 Suj„ X -t- 60 yf x = 0 


(1-62) 



This is a linear differential equation with constant 
coefficients. Its characteristic equation is clearly 



+ + gu/ =o 



(1-63) 



Three cases arise in obtaining solutions to the 
differential Equation 1-62. 

As an example, we will take 

Xo « 4 ft (1-64) 

Xo = 32 ft/sec (1-65) 

Cass 1 : C r.it ic ally-Pant p . s d-MQ t ibn.- ( 5 _t 1 ). 

As an example we will take 

c = 9.52 Ibf /-ft. sec" 1 (1-66) 

so that 



5 = 3 . 



2m C On 



= i 



(1-67) 



In this case, the characteristic equation has a double 



root 



51 



= - OJy 



(1-68) 
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and the most general solution is 



x (i ) = (A+ + A.-t) e"^ (1-69) 

where A + and A_ are constants determined by the initial 
conditions as follows 

Av = (1-70) 



and 



A_ - x„ + X o 0tJ, 



(1-71) 



so that the solution in terms of the initial conditions can 
be written 

x(l) = [x 0 (if U) n i) + Ki] (1-72) 

For our example 

xOQ = [ 4 ( If 8.51 -t) + 32 i] o' 8 ''’** (1-73) 

or 

„ . , , - 8 .SU 

xd) - (4 + 66.041) e (1-74) 

Case 2; QY.gxdamp.ad_filQ..tlQJi Jllll 

As an example we will take 

c = 14.28 lbf/ft-sec ' 1 (1-75) 

so that 

>* c 

6 = 2m u)m ' 5 (1-76) 

In this case the roots of the characteristic equation are 
real and negative 

^ - %u)„ ± con C5 2 -i) lz 

and the most general solution is 



(1-77) 



where A + and A_ are determined by the initial conditions as 
follows 





x„ = A+ + A_ 


' (1-79) 




(A + + A_) fA*-A_) 


(1-80) 


which yields 






* 


Xo 

(/5M *5)x„* 




/-> 4. 


2 /5M 


(1-81) 


A = 


C/lM' -5 )x.-3; 






2/5? T 


(1-82) 


For our example 






A + = 6.36 ft 


(1-83) 


A _ = -2. 


36 ft 


(1-84) 


and 

rn -I2.77i 9.5 H -9.5Ua 

xU) = <2 (6.36 <2 - 2.56 <2 ) 


(1-85) 


or 

*te) = 


/ 

- 3.2Si -22.28 i 

6.56 <2 - 2.36 e 


(1-86) 


Case 3 1 Underdamped Motion (5<1) 




As an example we 


will take 




c * 4.76 


lbf/ft-sec' 1 


(1-87) 


so that 

5 - 


q = 0.5 

2rr)iVr) 


(1-88) 



In this case the roots of the characteristic equation are 
imaginary with negative real parts. 



Si 



± j C Or, /4-$ Z 



(1-89) 
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If we define 



tt^ . co„ A-S* C for 5 < 1 J 

the so called damped frequency of oscillations, 

S\ s ' 5u)* * j 
and the most general solution is 



-font 






-fodt 



xtt) = e C 3 ” + A. g / 

Expanding the exponential terms we obtain 



where 



and 



x(t) = e ( a cos a) d i sin a) d t ) 

<* A+ + A. 

/S - j f A* -A_; 



The solution can be formulated as 



/, » \r t 

*(t) = X e cos C coj t f 0„ ) 

where X and Q a are constants determined by the 
conditions as follows: From Equation 1-69 we obtain 

e^°” x(t) =• X cos Ccu^i+Go) 

Differentiating Equation 1-97 with respect to time we 



For 



€ n [ 5u>„ x(t) + x(t)] = -X (v d smCoytfGo) 

t=0, Equations 1-97 and 1-98 become respectively 

x 0 =X cos ©o 



(1-90) 

(1-91) 

(1-92) 

(1-93) 

(1-94) 

(1-95) 

(1-96) 
ini ti al 

(1-97) 

obtain 

(1-98) 

(1-99) 
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and 



( 1 - 100 ) 



x„ + x 0 = -X ay sm 0 O 

So the constants are determined as 



0o - Ian 1 I 5ii)r) x 0 +• — — — 1 
*- — (Oj x. -* 

and 



( 1 - 101 ) 



x .i[ 



2 

X D t 



( ^UJyiXp <- X„ Y 



i 1/2 



( 1 - 102 ) 



Note: In Equation 1-101, keep the sign, so that the 

quantity in brackets will give the right quadrant for the 
angle & . 

For our example 

u) d = 8.5( Oo.5*)' /2 =7.37 sec' 1 (1-103) 

0O = 121.02° *> 2.1! rad (1-104) 

X » ± 7.76 ft (1-105) 



so that the response is 

/.i -4.261 

x(t) =-7.76 <2 COS (7.371 + 2.11 rad) (1-106) 

The roots of the characteristic equation can be plotted 
in the complex plane as shown in Figure 1-5. 

In cases 1 and 2 f the damping is so large that no 
oscillations take place and the mass simply returns gradually 
to the equilibrium position x=0. In those cases the roots 
are real and negative. 

In case 3, the damping has been reduced, so that 
oscillations about the equilibrium position, do take place, 
the amplitude of which dies out as t -*co . 
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* 



Hi 



•Figure 1-5 

Plot of the Characteristic Roots in the Complex Plane 



Figure 1-6 shows the responses in the above three 
described cases where it is assumed that the initial 
displacement is , and the initial velocity is x»=0. 




Figure 1-6 

Response of Displacement (x) of a Spring Mass Damper System 
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Notice that in the case where the damping 5 is zero, we 
have a spring mass system, without friction, described in 
section 1.4. In this case, the roots are imaginary and is 
the same as w-n (Figure 1-5) . 

It is clear that the system's physical constants (in our 
case 5,m,c,k), uniquely define the roots of the 
characteristic equation. Those roots are given the name: 
normal modes because one can determine from those roots, the 
corresponding characteristic motion (in our case, undamped 
critically damped, overdaraped, underdamped) . 

We have seen so far how an harmonic damped oscillator 
responds, if initially disturbed from its equillibrium 
position. Then knowing the initial conditions we can always 
find the response of the displacement or velocity versus 
time . 

Now we are interested in what happens if the oscillator 
is driven by an external input force F(t). 

In this case the solution will consist of two parts: 

xll) * x hom tt) ^ Cl) (1-107) 

where *vioy»>W is the homogenous part of the solution 
discussed so far, i.e. considering the input force zero and 
y pur-l (1) is the particular solution due to this particular 
applied force, i.e. considering the initial conditions zero. 

To obtain particular solutions one can use the method of 
undetermined coefficients, the method of variation of 
parameters or the Laplace Transform method. 
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In illustrating the above methods we are considering the 
following cases: 

(a) Sinusoidal Inputs 

The spring mass damper is driven bythe sinusoidal force 
F ( t) = F 0 cos co t (1-108) 

As an example, we will take 

F(t) = 24 cos 3.5 t (1-109) 

The method of undetermined coefficients assumes a 
superposition of sinusoidal solutions of the form 

*part W = A + coscot +• A_ sm cut (1-110) 

or equivalently 

x port ^ = X cos ((jut+ 0 O ) (1-111) 

where X and 0. are the coeffieients to be evaluated. 

Substituting into the differential equation 

x + 25to n x + o>yf x = F 0 coscot (1-112) 

we obtain 

C tu-n - oo* ) cosCuut+Oo) - 25tOnto sin (tot *-0<>) = ~ cos tot (1-113) 

Expanding cosine and sine and collecting terms, we get 

[ (ton -to 2 ) cos 0 O - 2 $u> n (o sin 0 O ] cos to t 



-[(U>n 2 -C0*) Sin 0c. + 25tOntO COS0o] smcot = ~ cosiot 


(1-114) 


Equating coefficients 




p. 

(air? -to*) cos 0 O - 25 u ) n (0 sm 0„ 

-X 


(1-115) 


and 




C o>y, 2 - oo 2 ) sm 0 O + 2$u>r,co cos 0 O = 0 


(1-116) 
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From Equation 1-116 we obtain 



. n -25 OOntJO 

-tan U 0 = — — — -7- 

Wrf -CO 2 



(1-117) 

Multiplying Equation 1-115 by cosSj, Equation 1-116 by 
sin f and by adding, we obtain 



C COy) 2 - CO 1 ) = -Z- COS 0 O 



(1-118) 



The cos 9c is expressed in terms of tan as 
COS 0 O = 



Vian0® tOn-U) 2 






(1-119) 



so that 



and 






F„ 



L(^-U)^ 2 t(25(0„co) 2 ] </2 



0 O - icin' 1 



-25oOrTCO 

0J fl 2 -C0 2 



( 1 - 120 ) 

( 1 - 121 ) 



For our example with the three different 
considered, the responses are: 

-for $ r i XpciriC-O = 0.28 cos (3.5 { +5.5 rad) 
for y pa,! (t) = 0.22 cos ( 3.51 + 5.21 rad ) 

-for 15~0.5 y par{(i) =0.56 cos (3.51+5.82 rcid) 

The total responses correspondingly are: 



values of 5 



( 1 - 122 ) 

(1-123) 

(1-124) 






x(U = 



-8.51 1 



(4+6&.041 ) <2 

+ 0.28 cos (3.51 + 5.5 rad) 



. z - -3.261 - 22.28 i 

6.56 e - 2 - 36 q 



+ 0.22 cos ( 3.5 t +5.34 rad) 



(1-125) 

(1-126) 
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(1-127) 



, , -4.26 i 

x[i) = -7,76 e cos (7.57 -t + 2.H racl) 

4 0.56 cos C 5.5 -t +5.32 rad) 

(b) Singularity Inputs 

The Laplace Transform method can be very easily applied, 
since we are dealing with linear differential equations 
having constant coefficients. 

Appendix A gives a short Laplace Transform Table and some 



Laplace Transform Theorems that help 
differential equations to algebraic. 


to transform 


In illustrating the above method, we are considering unit 
impulse and unit step input forces. 

(1) Unit Impulse Response 


The spring mass damper is driven by 

input 


the singularity 


FCO = w 

The differential equation becomes 


(1-128) 


x 4 25uj n x t a>n 2 x = S(t) 


(1-129) 


Neglecting initial conditions 
transforming, we obtain: 


and Laplace 


Xpari ^ = ( S 2 + 25co n s + U)n ) = 1 

or 


(1-130) 


Xpar-l (^) = ^ / S 2 +• S +■ cOri 

For our example with 5=1 


(1-131) 


Xparl( s ) = VCs + 85i) 2 


(1-132) 


, -8.51 -t 

x par-l CD = i e 


(1-133) 
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(1-134) 



and the total response is 



x(t) = (4 + 66.04 t ) e 8 ‘ 5,i + i g' 8 ' 5 ' i 



with 5 =1 .5 



Xpert (s) = Vs 2 + 25.55s <-7 2.42 = 4 /Cs+3.26)Cs+ 22.28) = °^ 5 iL y °: 052 . 5 - 

s.5.26 SV22.28 (1 _ 135) 

Xpart Ct) * 0.0525 e" 326i >0.0525 e ' 22 ' 28 ^ (1-136) 

and the total response is 

y(0 = 6-56 <f i2Qi - 2.56 e 2Z ' 2Si +0.0525 e 3-76 1 - 0.0525 e 22 ’ 28 * (1-137) 

With 5=0.5 



Xpert C s ) = V S 2 + 8.51 s +¥2.42 = VCs+ 4.26 +j ?.$7)( s + 4.26 -j 7.37) 



0.07 e 



j90’ 



0.07 <? 



-i 90° 



5+ 4.26 +j 7.37 5+ 4.26 -j 737 

x(i) = 0.14 e 42Gi cos (7.571 + t.57 rad) 



(1-138) 

(1-139) 



and the total response is 



x(i) =-776 e ' cos (7.37 1 + 2.U rad) + 

-4 

+ 0,14 e cos (7.37+ t.57 rad) 



(1-140) 



(2) Unit Step Response 

The spring mass damper is driven by the singularity 



input 



F(U -- u (£) 

The differential equation becomes: 

X + 2 ScOrr x + (0„ 2 X = U (1) 



(1-141) 



(1-142) 

Neglecting initial conditions and Laplace transforming, we 
obtain : 

Xpert (s) - U 2 + 25uj n s+a)n 2 ) = s (1-143) 
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or 



Xpart (s) = 5 /s 2Sa)r,s + co n 2 
For our example with 5=1 



8.51 



X pari ( s ) - Cs-t-8.51) 2 = s^8.S< " Cs+8.50 2 

-8.SK -8.51 1 

Xpari(i)*e - 8 . 5 H e 



and the total response is 

x % -8.51 1 -8.51 i - -8.SH 

x(l) = (4 + 66.041) e + e • i s 



With 5=1.5 
Xpcirl ( s ) = 



0.17 



-4.1? 



( s+S 26 ) Cs+ 22 . 28 ) S+ 5-26 S+ 22.28 



Xpart tO = 0-17 e 



• 3.26t 



- U? <2 



- 22.281 



and the total response is 
-3.261 , 

= 6.56 c ~ 4 

With $ =0 .5 



,.x ... -5.261 a -22.281 -3.261 -22.231 

x(t) = 6.56e -2.56e +0.17 e -U?e 



0.58 e 



-30.11° 



(1-144) 



0.s8e 



30.H e 



(1-145) 

(1-146) 

(1-147) 



(1-148) 

(1-149) 

(1-150) 



/pcirtO) - 1.57 ) (S+ 4 . 26 -j 157 ) ' 5 + 4.26 + j 7.57 * 5 + 4 . 26 -j 7 - 3 ? ( 1 - 151 ) 



Xpari^) - 1 46<3 426i COS ( 7.57 1 <• 0.55 rad) 



(1-152) 



and the total response is 



x(l)= - 7.76 2 4 26 1 cos (7.571 + 2. U rad) + 



-4.26 1 

+ 1.16 e cos (7.57 1 +0.55 fad) 



(1-153) 
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1.6 AXES SYSTEMS 

In describing the motion of our flight vehicle, we will 
make the following two assumptions: 

First, we will establish an inertial frame of reference 
on the earth considering the earth to be fixed in space. In 
this inertial frame, Newton's laws are valid. 

Second, the vehicle will be assumed to behave like a 
rigid body so that translational motion can be described by 
considering it as a particle with all its mass located at the 
center of mass, and rotational motion can be described by 
considering moments about that center of mass and constant 
moments of inertia. 

To describe the translational and the rotational motion, 
one has to establish a second non-inertial axes system, fixed 
on the vehicle with origin at its center of mass. In this 
reference frame Newton's laws do not hold, in general, and we 
have to add additional terms to our basic fundamental 
equations. 

It is convenient to use cartesian coordinate systems for 
both reference frames used. 

Let's consider the translational motion first. We will 
assume that 0xy2 is the inertial reference frame and 0 x'ijV 
the non-inertial. As in Figure 1-7 shown, r and r are 
position vectors pointing to an object and originating from 
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the origins of the coordinate axes 
denotes the position vector of the 
measured from the inertial origin. 



correspondingly. r 0 
non-inertial origin 




Figure 1-7 

Inertial and Non-inertial Reference Frames 
Translational Motion 

The relation between those three vectors is clearly 

r r r + r 0 

By differentiating twice, with respect to time, we obtain 

7 = r' 4 r 0 (1-154) 

Now since we have assumed that Newton's laws hold on the 
frame Ox yz , we may write 

m r = F (1-155) 

or substituting Equation 1-154 

W 7 r F - yr)7 e (1-156) 

where on is the mass of the object and F the force applied 
onto it. 
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Equation 1-156 merely says that we need the extra term 

"mr 0 to describe the translational motion of the object in the 

non-inertial frame. 

* * r / 

If r 0 = 0 then the frame OxtjV is inertial and Newton's laws 
hold. In this case. Equation 1-156 takes our familiar form 

nrt ?' = F (1-157) 

which is similar to Equation 1-155. 

Let's now consider the rotational case. We will assume 
that r o =0 so that o' is the same as 0 , as shown in 
Figure 1-8. 




Figure 1-8 

Inertial and Non-inertial Reference Frames 
Rotational Motion 



The position vector r can be expressed in the inertial 
frame as 



A A . 

rsxi+yj+zk 

or in the non-inertial frame 

A 

/ A /A / 1 

r — x i4.yj+-z k 



(1-158) 

(1-159) 
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where (*,y,z) and ( x' , y' , z' ) are the corresponding axes 
components of r . 

Similarly the vector A can be expressed as 

A = A* t + Ay £ + A z k (1-160) 

or 

A ^ Ax f Ay j' + Az k' (1-161) 

Differentiating Equation 1-160, 1-161 with respect to 
time, we obtain 



dA dAx A d Ay * dAz * 

di = d T 1 + "dT J ^ ‘dT k (1-162) 

and /dA\' dAx w dAij »/ cJAz a/ 

l dW " di 1 + cH * + * (1-16 3) 

where ( dA/^ )' denotes that we are taking the time derivative 

of A on the primed non-inertial system where the unit vectors 

A / A / A / 

i , j , k are constants. 

We are interested now in the time derivative of the 
vector A relative to the inertial system, but we would like 
to differentiate Equation 1-161 instead of Equation 1-160, 
obtaining thus: 



dA dAx a' dAi; a' dAz f 

dT"- IT 1 + dT J + dT k 




1 ' dr 

dT + 



a' dk 
Az dT 



or 




(1-164) 

(1-165) 



It can be shown that a more convenient expression for the 
three last terms is given by 
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= <jj * A 



(1-166) 




r 



dt 



+ 



' dj 

* dT 



a; 




oH 



where u) is the angular velocity vector of the primed system 
which is assumed to rotate about an axis having the direction 
of <£> , through the origin, and with angular velocity u>, 
relative to the unprimed inertial frame. 

For proof, assume for a moment that ^ is a vector "at 
rest" on the primed system so that 

=0 (1-167) 



We shall show that 

d Vdfc = w x A 



(1-168) 



Using the time derivative definition, we can write: 



dA p AU4-A t)-A(U - r 

dri: 0 — m — d-169) 

Figure 1-9 shows two successive positions of the vector A 
which rotates with angular velocity oo with respect to the 
inertial frame, so that the tip of the vector traces a circle 
perpendicular to the direction of w. 

To show Equation 1-169, we must show that the directions 
and magnitudes of the vectors on the left and right-hand side 
are the same. 

Notice that the- direction of the vector olA/c^ , by 
Equation 1-169, is the same as the direction of the vector 
A(-(.+ At)- Ml) , which is the same as the direction of the 
vector u>vA by looking at Figure 1-9. 
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Figure 1-9 

Time Derivative of a Rotating Vector 



The magnitude of the f!e£t-hand side of Equation 1-169 

is: 

I A (1+ A-t) - A (01 (AsinttOCwAi) 

5 mZ Ti (1-170) 

or 



.im 



dA 

di 



A (xj sin<fc> 



(1-171) 



which is the same as the magnitude of the right-hand side 
vector , I w y A I . 

Hence, both directions and magnitudes are the same, so 
Equation 1-168 holds and if A is not stationary with respect 
to the non-inertial system, we have according to 
Equation 1-165 
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cM 

ctt 



/ 



(1-172) 




+ oo x A 



which says that the time derivative of a vector A relative to 
an inertial frame equals the time derivative of the same 
vector measured at the non-inertial frame plus the term o>*A 
where a> is the angular velocity of the non-inertial as seen 
by the inertial frame. 

Combining translational and rotational motions, one can 
easily arrive at the equation: 




(1-173) 



where r. is the linear velocity of the non-inertial with 
respect to the inertial frame. 

One step further, differentiating once again 
Equation 1-172, with respect to time, we obtain 




or 



d^A 




+ 2 w * 






(1-175) 



which says that the second time derivative of a vector A 
relative to an inertial frame equals the second time 
derivative of the same vector measured at the non-inertial 
frame, plus some extra terms involving , A , du>/dl , and 
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(dA/di ) where u> is again the angular velocity of the 
non- inertial as seen by the inertial frame. 

Combining translational and rotational motions one can 
easily arrive at the equation 







(1-176) 



where r 0 is the linear acceleration of the non-inertial with 
respect to the inertial frame. 

A special significance of Equations 1-173 and 1-176 is 
given by substituting A with the position vector f of an 
object. They are taking now the form 



dr 

dT 



dr - - ~ 

— — + u)xr + r„ 



or 



V = (V)' + ujxr + r; 

— * 

where V is the velocity of the object, anc 



(1-177) 



(1-178) 




fo 



(1-179) 



or 



a = Ca )' -t- 2 <3 x ( V ) 4- uj x ( u)xr ) + dSxr + d„ (1-180) 

where a is the acceleration of the object. 

One can recognize that 

2 (3 x C V Y = corriolis acceleration (1-181) 

and exists if the object has a velocity (V) 7 with respect to 
the non-inertial frame. Also, 

uix(ujxr) = centripedal acceleration (1-182) 
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which is always directed toward the axis of rotation and is 
perpendicular to that axis, as shown in Figure 1-10. 




Figure 1-10 

Centripetal Acceleration 

Since we have assumed that the unprimed system is 
inertial, we may write 

oma = F (1-183) 

and substituting Equation 1-180 

omfd)' - F - 2in[(3x('/) , J -if?[$x(u>vr)J-ori(ujxr) - r ma 0 (1-184) 

where is the mass of the object and F the force applied to 
it. 

Equation 1-184 merely says that we need four extra terms 
to describe the translational and the rotational motion in a 
non- inertial frame. 

It is clear that the forces: 



46 



(1-185) 



2-m [ tu x ( = corriolis force 
om [ dix ( cox r)J = centripedal force (1-186) 

'Yn(uixr) = Euler force (1-187) 

<inna» = D 1 Alembert force (1-188) 

are not true forces. They only exist due to the fact that 
our object is moving on a non-inertial frame. 

According to our first assumption that the earth is the 
inertial frame of reference , we consider that it is rotating 
with fixed angular velocity a> so that we can write 
Equation 1-184 as: 

'rn(a)' = F ~ Zm a) x (V) - oox(u>xr ) - W cu (1-189) 

We have seen, up to now, the effects of considering a 
non-inertial frame of reference and our next step will be to 
orient our flight vehicle in that frame. The reason for 
selecting a non-inertial frame will be greatly appreciated 
when one will see how easy it is to formulate equations of 
motion of a body attached in this frame and moving with it. 

For flight control applications, it is appropriate to 
select the origin of the axis system on the center of mass of 
the vehicle. 

From then on we will orient the nose of the vehicle 
towards the positive x-axis, which will also be a symmetry 
axis. The positive y-axis will point to the right and the 
positive z-axis will point downwards having the direction of 
gravity, if one considers straight level flight as 
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Figure 1-11 shows, so that the xz plane will be a plane of 
symmetry. 

Vector quantities associated with this axis system are 
shown in Figure 1-11 and are described as follows: 



VERTICAL 
AXIS 



M 

% 

r 



positive pitching 
pitching velocity 

aerodynamic -force 
velocity 




L positive rolling moment 
p rolling velocity 

X aerodynamic -force 
o velocity 



M positive yawing moment 
r yawing velocity 
Z aerodynamic force 
w velocity 



Figure 1-11 
Axes System 

Velocity of the center of mass with respect to the 
inertial frame measured on this axis system 

V rut-f-vj+wi (1-190) 

Angular velocity about the center of mass 

cu = pi + r t (1-191) 
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(1-192) 



Aerodynamic force acting on the center of mass 

? -XI + Y j ,Zk 

Moment about the center of mass 

G = Lt + M? fNk (1-193) 

Directions as shown in Figure 1-11 
Also, we can define: (see Figure 1-12) 

Angle of attack: 

< x = -lari’ 1 ( w /o) (1-194) 




Angle of sideslip 



jS = sitT' (Vv) 



(1-195) 



The reference axis system has been adopted by B. Etkin, 
D. McRuer and others, and is usually referred to as body axis 
system (see references). The notation used is quite common 
on most references. 
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Another axes system adopted by C. Perkins is shown in 
Figure 1-13. It is usually referred to as wind or stability 
axes system. 




Figure 1-13 
Wind Axes System 

In this system the x-axis points to the opposite 
direction to that of the relative wind. Notice that in the 
undisturbed flight the body axes coincided with the wind 
axes, but in the disturbed flight that may or may not happen. 

1.7 RIGID BODY DYNAMICS 

Rigid body is a system of particles in which the distance 
between any two of them is not changing regardless of the 
acting forces. For our flight vehicle the rigid body 
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assumption implies that as far as its motion is concerned in 
three dimentions, only negligible errors are introduced. 

We will consider first the planar motion of rigid body 
where we will find the kinetic energy and the angular 
momentum expressions. 

Assume that the rigid body in Figure 1-14 can move in the 
xy-plane and can rotate about an axis perpendicular to that 
plane passing through the center of mass c. 




Figure 1-14 

Rigid Body - Planar Motion 

We denote ^ , the position vector of the center of 
mass c, 7 the position vector of any particle on the rigid 
body and r> the vector from the center of mass to that 
particle. 

Usually, the mass of the rigid body is not uniformally 
distributed over the entire volume V. It is convenient to 
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think for a moment that the entire volume of the rigid body 
containing a mass on, is projected on the xy-plane on an 
area A. Let's denote the mass per unit projected area of 
the body. It is clearly seen that the summation over the 
entire area of all those masses ^dA gives the total mass of 
the body. This mathematically is expressed as: 



orri = J 

From Figure 1-14 we have: 



c\A 



r =• r 0 + r 



(1-196) 

(1-197) 



For a differential mass ^dA at which the r vector points, 
the kinetic energy is expressed as ^ r 2 . integration over 
the entire area A gives the total kinetic energy of the rigid 
body. 






■ij>" 



(1-198) 



The quantity f 2 according to Equation 1-197 is expressed 



as 



r 2 = r-r - Cr; + r')-(r Q + r') = r„ 2 + r 2 +2r„r' 
Then Equation 1-198 becomes 



(1-199) 




The last term is zero since 




= r 0 yn-r a m *o 



( 1 - 200 ) 



( 1 - 201 ) 
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We now express the velocity f' as a sum of two 
velocities : 



r 




+ r 0 9 



( 1 - 202 ) 



where dr'/dl is the radial velocity, 0 is the angular 
velocity of the body, r'0 the tangential velocity of the 
particle , r and x) denote the corresponding unit vectors in 
polar coordinates. 

Noticing that dr'/c)l=0, Equation 1-200 is written 



We define 




(1-203) 



I = f r^dA 

•* A 



(1-204) 

the moment of inertia of the body about an axis perpendicular 
to the plane of motion through the center of mass. Then 



K = 2 + 2 I 0 (1-205) 

This equation merely says that, in addition to the term 
| • ^ 

- mf, that describes the translational motion, we need the 

z 

extra term -IQ 2 which describes the rotational motion of 
2 

the body around an axis perpendicular to the plane of motion 
through the center of mass. 

The angular momentum of the same differential mass vjcJA 
about the origin is expressed as 



r x 



*1 



dA 



(1-206) 
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Integration over the entire A gives the total angular 
momentum of the rigid body. 



h - j C r x r ) dA 



The cross product term can be written 



(1-207) 



r vr = (£ + r')x (£ + r') = r 0 x f; + rxr' t r'x r„ +r 0 xr' 
Then Equation 1-207 becomes 



(1-208) 



k = r.xr; | n dA + j (rVr'J^dA + r a x r'rjdA * r 0 x J r'rjdA (1-209) 

In this expression the integral of the last term is zero 
as Equation 1-201 shows. By a similar argument the third 

. * A 

term is also zero. Noticing that r / =r0? , Equation 1-209 
becomes 

k =(r 0 y r;)w + 0 j r' ( r'x (1-210) 

— * 

If we use the center of mass at the origin ( r . =0) we are 
left with: 



k = 0 f r(r x§) (1-211) 

sA 

The cross product in the integral indicates that the 
angular momentum vector is perpendicular to the plane of 
motion. Since I r x D ( =r f the magnitude of the angular 
momentum is given by: 



k 




r^d A *10 



( 1 - 212 ) 
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So the magnitude of the angular momentum of a rigid body 
rotating in a plane, about an axis, through the center of 
mass, is given by Equation 1-212. The angular momentum 
vector points in the direction of the axis of rotation that 
can be easily found by the right hand rule. 

In summary, we have seen that the motion of a rigid body 
in planar motion can be described by a translational motion 
of the center of mass, plus a rotational motion about an axis 
through the center of mass. 

We will now look at the motion of a rigid body in a three 
dimensional space. It generally has six degrees of freedom 
if not constrained. Hence, we need six coordinates to 
specify the motion. The three of them can describe the 
translational motion of the center of mass of the body and 
the remaining three, the rotational motion. 

We consider two coordinate systems. An unprimed system, 
as shown in Figure 1-15, whose origin is 0, for simplicity is 
fixed on the rigid body but its axes do not rotate with 
respect to an inertial frame, and a primed system whose 
origin and axes are fixed on the rigid body. 

We also assume that the rigid body has a mass 
distribution per unit volume (density) which, in general, is 
a function of the position vector f to a particle on the 
rigid body having a differential volume cV . The total mass 
is given by: 
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Figure 1-15 

Rigid Body - Space Rotational Motion 






-I p f?,d ' 



(1-213) 

where p(r) is the density as a function of r and dy a small 
differential volume where r points to. 

For the differential volume dx the angular momentum for 
rotation about the origin of the unprimed system is given by 

r x [ p(r)dv -F] (1-214) 

Integration over the entire volume gives the total angular 
momentum: _ 

k =!_(?*?) P (?)dv (1-215) 

where f' , r are measured on the unprimed system. 

Let oo the angular velocity of the rigid body about the 
origin of the unprimed system and r the position vector from 
the primed system. 
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It is more convenient to express F , which can be seen as 
the velocity measured on an inertial frame, as 

F = F' + aJ x F (1-216) 

Notice that 



/ 

r 



dr' ^ , 



r'0 V 



(1-217) 



where £' and t) are unit vectors in polar coordinates of the 
primed system. 

It is clear that dr*' /<£ =0 according to rigid body 
assumption, and §=0 because there is no relative rotation of 
the differential volume civ that can be seen by the primed 
system. In other words, F' is a vector at rest in the primed 
system so that we can write 



# 

F = 65 x r 


(1-218) 


Substituting in Equation 1-215) , we obtain 
h =• j* r x ( ujx r ) p(r) dv 


(1-219) 


Using now the tripple product identity 
A * (sxc) = (A-C)i - (A -B)c 


(1-220) 


we can write 

r x (Jxr) - - ( cF • r ) F 


(1-221) 


so that Equation 1-215 becomes 

It = | £ r z - fo3-F) r J p(r) dv 


(1-222) 


In component form we can write 

h L = j* [ f z — (^ d}- r ) ] p C?) civ 


(1-223) 



where 1=1, 2, or 3 denotes the x, y, or z component 

respectively. The dot product is expressed as: 
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(1-224) 



CD 



? = u;,r 4 * u> 2 r, + c o 3 r 3 r j( Dj 



That permits us to write 

K 






(1-225) 



(1-226) 



i = L \ [ r *&'j - ^ r j ] Wj ptr) civ 

J--1 K 

where cij =1 if i=j and ^ij =0 if i£j because 

3 

/ . ^'lj = ^ij +■ ^12 t^ 2 f S15 tDj = CD • 

J--< 

We can now define 

J 'i =i v ( ri 5ij -r i0 ] pfrjW^ (1 . 227) 

as the ij th component of the inertia tensor I about the 
origin of the unprimed system. 

So finally, we express the i th component of the angular 
momentum as : 



3 

r — _ 

^ = Z- ^ W j 

J=1 



(1-228) 



In the expanded matrix form this gives: 



k 




S 


II 


k 





I 





ky 


I x r 




Iyy 


J y 2 


X2 


Jyz 


1 22 



uJ x 

(O y 
<V Z 



(1-229) 



This result is compared with the result of the planar 
motion case that is expressed by Equation 1-212, where the 
inertia scalar term is expressed now as a matrix. 

The components of the inertia tensor are given by 
Equation 1-227 in cartesian coordinates as: 
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(1-230) 



^ ' fff [ (Vh^+z*) -X 2 ] pCr) dxdijdz: = 

~fjj pC r ) dxdydz 

Similarly , 




( x 2 4 z z ) p(r ) cJ x. d y d z 



(1-231) 



zz 



*JJJ C x*+y 2 ) p(rj dxdydz 



(1-232) 



Also, 



JxlJ = 


pc?) dxdydz 


(1-233) 


Jxz * Izx XZ 


p(r) dxdydz 


(1-234) 



Example: Calculate the inertia tensor of a cube of 

side a, mass m and uniform density p about one of its corners 
(Figure 1-16) . 




Figure 1-16 

Inertia Tensor-Calculation of a Cube 
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J*x p[f( dxcJydz =r TL I I I (yt+L 2 ) cizdydx 

^ ^ 4? 'O -k 



Oftt 



a ,.ct 




)«M* ' r.ffT«‘ r,) dx . 



o -fc 



a s I ' 5 3 

J o 



a J V 3 



a \ ™ 2a 2 2 

T a ) = — -- n»a 



a 5 3 



(l 



Similarly, 



1 33 c Izz = T ' yrlCl2 



(!■ 



5 rfr,*v<- 



I«, r -p hy 



o J o J o 



<ta f tt f a I , nn f a a ? i 

* JiJ. . *'?J. *2 ady = 



m\ a J a 2 < 2 

— mi a 

a 3 Z 2 4 



Similarly, 



I xz 



— mo. 



(1 



(1 



so that the inertia tensor is written: 



I r tmo‘ 



2/3 -'/4 -% 

-*/4 2 /3 - V 4 

-</4 -74 2 /3 



(1 
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235) 

236) 



237) 
23 8) 

-239) 



The eigenvalues and the eigenvectors of an inertia tensor 
have a very special significance. The eigenvectors determine 
the directions of the so-called principal axes of a rigid 
body that pass through the point where the inertia tensor is 
calculated (usually the center of mass). The corresponding 
eigenvalues are called the principal moments of inertia. 

It can be shown that if a rigid body rotates about one of 
its principal axes, its angular momentum about the point 
which the inertia tensor is calculated, is parallel to its 
angular velocity. 

For proof, let $ a principal axis, and B the 
corresponding principal moment of inertia. The eigenvalue 
problem equation which has to be satisfied is: 



Having rotation about a principal axis simply means that 
<Z has the direction of the principal axis which is expressed 
as: 



where C is a constant. In fact, if /$ is unit vector, C is 
the magnitude of the angular velocity to. 

Now the ith component of the angular momentum becomes 



If* « 



( 1 - 240 ) 



or 



iij h -*/*i 



( 1 - 241 ) 



U) = C /S 



( 1 - 242 ) 




( 1 - 243 ) 



and finally. 



hi * B to i 



( 1 - 244 ) 



61 



In vector form 



k = B uJ (1-245) 

which proves that k is parallel to a). 

It can also be shown that 

1. Any plane of symmetry of a rigid body is 
perpendicular to a principal axis. 

2. Any axis of symmetry of a rigid body is a 
principal axis. 

The density p may not be constant. For example, if the 
xy plane is a plane of symmetry, that simply means that 

p(-*>y,z) ~ (1-246) 

and if the z-axis is an axis of symmetry, that simply means 
that 

p(-x,-y,z) = p ( x, y , z ) (1-247) 

3. Since the inertia tensor of a rigid body is 
symmetric, (i.e. Jxy = Jyx etc.), the principal axes are always 
orthogonal to each other (Important result from Linear 
Algebra) . 

4. If two or three eigenvalues of a symmetric tensor 
are the same, there exist two or three orthogonal 
eigenvectors associated with that eigenvalue. 

Example: The inertia tensor of a rigid body has been 

calculated with the reference point at the center of mass and 
has been found to be 
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-JO 



I = 



15 0 

0 25 

-10 0 



0 

30 



Find the principal axes and the corresponding principal 
moments of inertia. 

Solution: 

Let jj denote the eigenvector (principal axis) and A the 
corresponding eigenvalue (principal moment of inertia) . 

The following eigenvalue equation has to be satisfied 



from which 





- T J 


= A I 






A-15 


0 


10 


H 

O 


0 


A-25 


0 




10 


0 


A- 30 



= 0 



(1-248) 



(1-249) 



Ol5)(A*25)(A-30) + JO [-10 025)] = O 
^ 025) [ 015)030) -100] s0 
=> (A- 25) 035) CA-JO) =0 

Hence, the principal moments of inertia are: 
A, = 10 A 2 = 25 A 5 = 55 

To find the eigenvectors: 

For Aj =10 Equation 1-24 9 becomes 



-5 0 

0 -15 
0 



■10 






= 0 



(1-250) 

(1-251) 

(1-252) 

(1-253) 



(1-254) 



-10 



63 



0 

-20 



which yields 



" 5 i)<x + J) <)t - 2X Z 

- 45 3 , y =0 =» *0 

“20 ^ 1Z -0 =* Ji 4x - 2* tz 

The eigenvector that satisfied Equation 1-254 is 



(1-255) 



Sii « 



2 

o 

4 



(1-256) 



It is more convenient to express the eigenvector as a 
unit vector 

' Vf5 






0 

V6 



(1-257) 



Similarly, we find that 
0 

i - 1 

0 



assosialed A 2 = 25 



(1-258) 



and 



Si, 




assosiaied wi-ffi A 3 - 5 5 



(1-259) 



It is important to note that the above eigenvectors are 
orthogonal since they are coming from a symmetric tensor. 

Since now the principal axes of a rigid body are 
orthogonal to one another, or can be chosen orthogonal, it is 
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very convenient to have them coincide with the primed system 
axes, the origin of which can be moved to the origin 0 of the 
the unpriraed system. The importance of this arrangement will 
be shown a little later. 

The kinetic energy of the differential volume dv is 
2 ( r • f ) p ( r ) dv. Integration over the entire volume gives the 
total kinetic energy: 




(1-260) 



Substituting 



r 




(1-261) 



we obtain 




(1-262) 



Using Equation 1-224, 1-226 and the vector 



identity 



( axb)(cxd) = (XcK8-5) - (A.5H8-C) 



(1-263) 



we can write 



(uixr }(u>xr) = r z 





(1-264) 



so that 




(1-265) 



and finally 



65 



(1-266) 
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Notice that if the body is rotating about a principal 
axis ( kill*)) its kinetic energy is maximum Ksjkoo. That 
happens when the rigid body rotates about a principal axis. 

This compares with the result of the planar motion case 
that is expressed by Equation 1-205, if we neglect the 
translational motion. 

We are now in a position to formulate the rigid body 
equations of motion. 

The translational motion equation is given by the 
equation of motion of the center of mass of the rigid body 

W C - F (1-267) 

or 



'vnf — + a>xV) =• F 

l olt ' (1-268) 

The rotational motion about the point (0) which is 
conveniently chosen as the center of mass (c) , as Figure 1-17 
shows, is given by 



dh 

eft 




(1-269) 



where G is the net external moment about the center of mass 
(c) and h is given by Equation 1-228. 
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z 



I 



i 

I 




Figure 1-17 

Rigid Body - Space Motion 



The vector dJi/dt is measured in the unprimed system, but 
it can be expressed as: 



dk 




(1-270) 



where ( dk/dl ) is the same quantity measured in the primed 
system, and tu is the angular velocity about the center of 
mass. 



So Equation 1-269 is written in component form as: 



(A) t (fixh)i - G, 

where i refers to a primed axis now. 
Equation 1-228 



(1-271) 

Or better, using 
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(1-272) 



where I lj is a component of the inertia tensor calculated 
along the primed axes. 

Notice that the primed axes unit vectors are vectors at 
rest, in the primed system, and obviously Iij is not changing 
with time. Hence, 



J-i 

This result is interpreted as follows: Suppose the rigid 

body rotates about a principal axis with constant angular 



moment required for this rotation is zero. 

This compares with the planar motion case result which 
comes by taking the time derivative of Equation 1-212. 



Equation 1-27 3 is the equation of rotational motion in 
the general case for rotation about a random axis, where 
external moment is required even for constant angular 
velocity. 

It is now clear why we formulate the equations of motion 
in the primed, non-inertial coordinate system. We have 
avoided the time differentiation of the inertia tensor in 
Equation 1-272 by paying a little for the extra term u)xl*. 




(1-273) 



velocity o3 . Rotation about a principal axis implies ui 
parallel to ^ , so that t&xla-g. This simply means that the 



k r Icb r G 



(1-274) 
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In summary, we found the equations of motion of the rigid 
body described in the convenient form which will be used to 
derive the flight vehicle's equations of motion in Chapter 2. 



They are 




(1-275) 



and 




(1-27 6) 



where 



k = Lj <Dj 



(1-277) 



and I ij is calculated about the center of mass with respect 
to the primed axes. 

Finally, if the primed axes are principal axes of the 
rigid body it is not difficult to see that the inertia tensor 
has only diagonal terms representing the principal moments of 
inertia (eigenvalues) . This result comes from the 
diagonalization process of a square matrix. So, in this case 
Equation 1-276 is simplified as: 




(1-278) 



where 



k - lit 



(1-279) 
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CHAPTER 2 



VEHICLE 1 S EQU ATIONS OF MOTION 
2.1 INTRODUCTION 

With the methods of obtaining rigid body equations of 
motion, established in the preceding chapter, we now proceed 
in this chapter to formulate the equations of motion in a 
most generally useful form. 

The equations are basically simple in the sense that they 
are derived by making use of the basic Newtonian laws for 
forces and moments. Great advantages are obtained, as we 
have seen, by using a non-inertial frame of reference, 
because with an axis system fixed in the body, the inertial 
properties are constant and that simplifies the moment 
equations. 

Some penalties are obtained by introducing non-inertial 
frames, the most important of which is the addition of extra 
terms in the basic equations of forces and moments, in trying 
to make the basic laws valid in this environment. 

In any case, although the equations are simple, they are 
highly non-linear and not generally solvable mathematically. 
The formulation is greatly simplified by the use of 
linearization techniques, and small pertubation theory. 
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This method assumes an equilibrium position which is, in 
fact, a solution to these equations. That allows for the 
changes in forces and moments to be expressed as linear 
functions, in terms of small pertubations in linear, angular 
velocities, and other parameters involved. 

The equations of motion will be derived for the body axis 
system established in Figure 1-11 and will be left in their 
diraentional form. Non-dimensional form will be mentioned as 
it applies to body axis system (Etkin's) and wind axis system 
(Perkin's) . 

Finally, the equations of motion will be split into two 
separate groups, the longitudinal and the lateral groups of 
equations. 

2.2 FORCE EQUATIONS 

Starting from the vector force equation 

IF +SyV) 

where 

IF = IX* + ZYj + 11 k 
V - 0 1 fVj-fwk 
i 

we can derive the force equations in component form. 

The cross product term is expanded as: 
ujx V = (cjw-vr) t +(-pw+ur)j+ k 




( 2 - 1 ) 

(2-2) 

(2-3) 

(2-4) 



(2-5) 
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( 2 - 6 ) 

(2-7) 

( 2 - 8 ) 



and the force equations of motion are simply 

2 X " = oft) ( b + - vr ) 

I Y - on ( v -pw -f-ur) 

IZ - -m (w t -pv -c^o ) 

These equations describe the translational motion of the 
center of mass. 

Consider the motion of the center of mass in the xz plane 
as Figure 2-1 shows, in two successive different locations on 
the flight path during time St. 





Motion of the Center of Mass in the xz-Plane 



The acceleration along the x-axis is expressed as: 
(im-6u)cosS 0 t -v 



CX x = 



It " (2-9) 

For small pitch angle 0, cos&9 =1 and sin60 =S0 , so that 

U+Su + -t-SwSG-u 



CU = 



St 



( 2 - 10 ) 
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letting now St^O and neglecting the higher order term <5 wjo© , 
we obtain 

a x =• v + (2-11) 

and for the corresponding force equation 

Z X = or) (u + wc^) (2-12) 

which is the same as Equation 2-6 by neglecting the term vr . 
Similarly, the acceleration along the z-axis is expressed 

as: 



a ; 



(w+ Sw) cosS9 - (u-hSu) sin 66 - w 

St 



which is simplified as: 



(2-13) 



a 2 



v^+Sw - uS0 - SoS© - w) 



St 



(2-14) 



or 



a z - w - cjv ( 2-15 ) 

and the corresponding force equation becomes: 

22 - onCw-cju) (2-16) 

which is the same as Equation 2-8 by neglecting the term p/ . 

If we now consider the motion of the center of mass in 
the xy plane as Figure 2-2 shows, in two successive different 
locations on the flight path during time St , we can derive 
the acceleration along the y-axis as: 



a 



i) 



(v+Svj cos Sip * (u+Sv) -v 

si 



(2-17) 
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Figure 2-2 

Motion of the Center of Mass in the xy-Plane 



For small yaw angle ip , cos5^=l and sin&4> = Sij>, so that 

(v+£v + -t-SuStp-v) 



a 



'i 



St 



(2-18) 



letting now <Sl-*o and neglecting the higher order term SuSy/ , 
we obtain 

a y =• v +vr (2-19) 

and for the corresponding force equation 

Z Y = on ( v +• ur ) (2-20) 

which is the same as Equation 2-7 by neglecting the term pw. 

In summary, the force equations are simplified neglecting 
higher order terms as: 



IX = 'vn ( V + Cjvj) 


(2-21) 


2Y = 'W ( y + tir ) 


(2-22) 


11 = oyt(w-^u) 


(2-23) 
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2.3 VEHICLE SYMMETRIES 



Up to this point we have made two assumptions. The first 
assumption was that earth is considered to be fixed in space 
so that we can establish an inertial reference frame on it. 
The second assumption was that the vehicle is considered to 
be a rigid body. 

Facts are simplified further by considering the 
symmetries of the vehicle. For an airplane or a missile, the 
xz-plane is a plane of symmetry, if one assumes that the 
following holds for the density distribution. 



Notice that Equation 2-24 says that p is an even function 
with respect to y but the integrant function is odd because 
it is a product of an even and odd function. Since the y 
integration has symmetric limits, the integral is zero. 

For the same reason 



The inertia tensor, therefore, for an airplane becomes: 



p (x, y,2 ) - p (Vy,z) 



(2-24) 



so that 




(2-25) 




(2-26) 



Jxx 0 -Ixz 

0 Iijij 0 

-Ixz 0 Izz 



(2-27) 
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Notice also that the y axis is an eigenvector of the 
inertia tensor (Equation 2-27) because it satisfies the 
eigenvalue matrix equation: 



Ixx 0 -Jy? 




0 




0 


0 Iy y 0 




J 


= hu 


i 


I XZ 0 Iz2 




0 




0 



Hence, the y-axis is a principal axis, with lyy being the 
corresponding eigenvalue or principal moment of inertia. 
This was expected since y-axis is perpendicular to the plane 
of symmetry xz. 

Unfortunately, this is the only symmetry that exists 
generally in any conventional aircraft. 

In addition to this symmetry another symmetry can be 
found in gun projectiles and in most missiles and 
space-crafts. They have a second plane of symmetry, the xy 
plane, so that one can assume 

p(x,y,z) = p (*»y,-2) (2-29) 

and hence, 

Ix 2 = I 2*«0 (2-30) 

Then, the inertia tensor for a two-symmetry vehicle 
becomes 



Ixx 0 0 

0 Iyy O 

0 0 J 22 



(2-31) 
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Notice that in addition to the y-axis being a principal 
axis, the z-axis, as well as the x-axis, are all principal 
axes with corresponding principal moments of inertia Izz , and 
Ix> , respectively. 



2.4 MOMENT EQUATIONS 



Starting 


from 


the vector moment equation: 








S -I f ♦ (5*iO 


(2-32) 


or in component form 




where 


Ci 


= 1 It, 3* *tSxCh 

J*l 


(2-33) 






= Iij Wj 


(2-34) 






a - Lt + Mj + N k 
k - hx t + hy j 4- h 2 k 


(2-35) 

(2-36) 






du) * A • . • | 

_ , f ‘ <-^1 * r t 


(2-37) 


Equation 


2-34 


is expanded as: 








k* =■ plxx - r Ixz 


(2-38) 






kvj « 

hz = -<p Ixz + r Izz 


(2-39) 

(2-40) 



also, the cross product terms is Equation 2-33 is expanded as 
( uj x k ) =(^h 2 -r(i y )l ^ (rh x --pk 2 ) J + (2-41) 
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or 



((ijxk) =(-|pI« t^ r l Z 2 - rc^ Iyy ) L 

■+ ( 7 Iyx - r 2 1„ + p 2 Ixz - rp hz ) j 
+ (^p I yy - <^p Ixx t Ixz ) k 



(2-42) 



or 



— ♦ 

( U) X It ) = £ + (jjf ( 1^2 - lyy ) J L 

4 [ rp (Ixx- 1 Z 2 ) + 1x2 (p 2 -H)J j 

+ [ p^ (Jyy-Ixx) 4 1 x 2 ] k 

Thus, the moment equations are: 

L = Ixx p - Ixz f - Ixz + Cjr (Iz2-Jyy) 
M - Iyy ^ + rp ( Ixx - I2Z ) + Ixz (p 2 -r 2 ) 

N - Iz2 r - Ixz p + p^ ( Iyy -Ixx ) 4 rc^ Ixz 



(2-43) 

(2-44) 

(2-45) 

(2-46) 



2.5 EXPANSION OF GRAVITY FORCES, THE LINEAR AND ANGULAR 
VELOCITY 

The total force acting on a flight vehicle can be split 
into a summation of aerodynamic forces (lift, drag and 
thrust) and a body force which is the gravity force. 
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Figure 2-3 shows an accidental orientation of the vehicle 
with respect to the vertical. The gravity force M can be 
expanded in components along the body axes as follows: 



X 

II 


- onj sin 0 


(2-47) 


w b - 


'witj cos 0 sin (J? 


(2-48) 




'Wq cos© cos 4> 


(2-49) 




Figure 2-3 

Expansion of the Gravity Force 

For the expansion of the linear velocity vector along the 
body axes consider in Figure 2-4 the inertial reference frame 
being y , 2^ . Transformation to the body axes coordinate 

system is obtrained by rotating the reference frame in an 
order-dependent manner as follows: 

(a) Rotate around reference axis 0z L by the yaw angle y>, 
which defines the intermediate reference system 0y 2 \} l z 2 = 

0*2 1*2 2 i • 
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(b) Rotate around intermediate axis 0 y 2 by the pitch 
angle 0 , which defines a second intermediate system 0x^,2^ = 
0xy 2 z 3 . 

(c) Rotate around final axis Ox by the roll angled, to 
define the desired body axes system 0xy2 = 0x 5 y2 . 




Figure 2-4 

Expansion of the Linear Velocity 

The components of the velocity vector in the body axes 
system (U,V,W) are defined in terms of the vector components 
in the fixed inertial frame ( U 4 ,V, ,W, ) by means of the angles 
and the intermediate vector components (Ui,V 2 ,W 2 ) and 
(U 3 ,Vj,W 3 ) . 
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The following relations hold: 
for the first rotation 

U 2 = U 1 cos ip + Vi sin ip (2-50) 

V 2 - -V l sinap + Vj cosip (2-51) 

= Wj (2-52) 

for the second rotation 

U 3 - 0 2 cos 0 - sm0 (2-53) 

V s = V 2 (2-54) 

U 2 sw0 + W 2 cos 0 (2-55) 

and for the third rotation 

U - U 3 (2-56) 

Y =■ V 3 cos4> 4 Wj sir) 4 ? ( 2 - 57 ) 



W = -V 3 sin 4? + Vl/j cos 4? (2-58) 

so that the velocity components along the body axes are 
expressed as 

U = Uj cosipcos© 4 sirup cos© - Wj si v6 (2-59) 

V - Vf (cosip sin© sin - sirup cos <p ) +y s (simp sir>0sm<i> -i- cosip cos 

4 W| cos0 simp (2-60) 

W = V S ( cos ip sin 0 cos <p + sin ip sw<p) + V, (simp sm© cos<fc> - cos tp sin 4> ) 

4 W, cos 0 cos (2-61) 

Similar approach can lead us from the body axes system to 
the inertial for which 

- V cos9 cosip + V (sm<P sw 0 cos ip -cos4>swzp) 

+ W(cos4> sm0 cosip 4 sm<p soup) (2-62) 



81 



V f = l) cos© simp + V ( smcp sin© simp + coscpcostp") 



+ W (cos<p sin0 sw jp - siricj? cos ip) 
W, = - V smG V V sir) <$> cos G -t- W cos 4> cos 0 



(2-64) 



(2-63) 



Integration of Equation 2-62, 2-63 and 2-64 give the 
expressions for the translational motion but in other than 
special simple cases they are not integrable. 

Following a similar procedure we can expand the angular 
velocities Prq f r into components along the body axes. 

With the vehicle in horizontal flight, one can see the 
yawing-velocity vector having the same direction as gravity, 
the rolling velocity vector having the direction of the 
longitudinal axis and the pitching velocity vector, 
perpendicular to both as shown in Figure 2-5. Those 
orientations are obtained according to right-hand rule. 




Figure 2-5 



Expansion of the Angular Velocity 
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Then the angular velocity vector cu can be expanded in the 
body axes as follows: 

V = - V Sin0 (2-65) 

© cos<f> -h ip cos© sm4? (2-66) 

I f = -0 sir>4> + %> cos© cos 4? (2-67) 

One can now solve the Equations 2-65, 2-66 , 2-67 for 0, 
V, $ . 

Multiplying Equation 2-66 by cos ij> , Equation 2-67 by 
sin t and subtracting, we obtain: 

9 r c^cos4> -rsin<J> (2-68) 

Multiplying Equation 2-66 by sivi4>. Equation 2-67 by 
cos <|? , and adding, we obtain 






sm<j> 



r 



cos 4? 



cos 9 cos© (2-69) 

Then substituting Equation 2-69 into Equation 2-65, we 
obtain 

4> = p + q sm4> -tan© * rcos^-tan© (2-70) 

Integration of Equations 2-68, 2-69 and 2-70 give the 
expressions for the rotational motion but in other than 
special simple cases they are not integrable. 

We can now modify the vehicle's equations of motion by 
considering the gravity as follows: 

Force equations: 

X = 'vn f i) + c^w -yr + <j sm© ) (2-71) 

Y - tr) ( v + vt - pw - cj cos 0 sm 4? ) ( 2-72 ) 

Z = /m ( + pv -yj cos © cos <p ) 
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(2-73) 



where now the X,Y,Z forces represent only aerodynamic forces 



(i.e., contributions of lift, drag and thrust). 
Moment equations: 


L = Ixx p - Ixz f “ C^p IxZ +• Cj_r ( l zz - lyy ) 


(2-74) 


M * lyy q + rp (lx* -lzz) + Ixz C p 2 -r 2 ) 


(2-75) 


N - lzz r - Ivz p + -pc^ (I^y-lxx) i-r^ Ixz 


(2-76) 


with auxiliary relations: 


'P - cj> - ipsm0 


(2-77) 


• • 

~ 9 costp +- ip cos© sirxt> 


(2-78) 


r - - 0 sir? 4> + cos ® cos ct> 


(2-79) 


Notice that the gravity force contributes 


zero moments 



because it acts at the center of mass. 

2.6 LINEARIZATION OF FLIGHT VEHICLE'S EQUATIONS 

The three force and the threee moment equations derived 
in Section 2.5 must be linearized for use in the stability 
and control analysis. It is assumed that the motion of the 
vehicle consists of small deviations about a reference steady 
flight path. 

A steady flight path is described by the general 
equations of motion, if we consider all time derivatives 
zero. We denote by capital letters and the subscript (o) , 
the steady flight conditions, obtaining thus: 

Xq = m - Vo£ 0 + sin 0o ) 

X o - wCUoRo - PoWo - g cos 0 O sin ) 

Z 0 - w CPoVo - QoV 0 ~QC os 0 O cos 4> 0 ) 
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(2-80) 

(2-81) 

(2-82) 



(2-83) 

(2-84) 

(2-85) 



Lo r ~ QoR) I>2 + C Izz- Iyy ) 

Mo - *oP, (Ixx-I M ) + Ixz (?*-&) 

No = 7oQ 0 C lyy -Ixx) + 1?«Qo 1*7 
By subtracting from the original equations the above 
steady flight equations, we obtain the perturbed equations. 
For example, subtracting Equation 2-80 from Equation 2-71 we 
get: 

X-Xo = w [ 0 + c^w - QoWo - vr +V 0 R 0 + (j ( sm0 - smGo)] (2-86) 

Notice that 

^ CQo+s^Hwo^w) -QoWo+w/o^ ♦ o a Sw j. ^r A0 ' r 

(2-87) 



Similarly 

VY = VoRo + Vo Sr + Ro 5v (2-88) 



Also 

sw0 - sinC0 o +60) - sin 0o cos 69 +cos0o 5inS0 (2-89) 

For small perturbation quantities. Equation 2-89 becomes: 

sin0 = sm0 o +(cos0 o )60 (2-90) 

Combining all these. Equation 2-86 becomes 

6X = /yn[ u-t- fQoSvu -V 0 6r -'PoSv +(^cos0»)3©J (2-91) 

Although we are deriving perturbed equations, it is more 
convenient to drop the differential symbol 6 . 

Similar work leads us then to the following set of 
linearized equations: 
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(2-92) 



X- W [ 0 +Wo^4-QoW - Vt> r - R 0 v <-(^jcos 9o) 6 J 
V = om [ V + U 0 ir + RoU - Wop - Vo w - (<j cos0 o cos <k) 4 > f(gsm6l£in^)9j(2-93) 

Z - on ^ w + Vo p 4- T ou - (Jo — Qo u +( cj cos 6o cos ) 4> -t-CgsDi&cosA*)®]^ - '^ ) 

L r Jxx p ~^xz r - (Qop +To<^)lxi + (Qor + Ro0 ? )(Izz *Iy</ ) ( 2-95 ) 

M« Iyy ^ + (Rop + ? 0 r) (Ivx-Izz) + 1« C 2P 0 p - 2fc> r) (2-96) 

M = J 2Z f - Ixz-p + (R»^4-Q 0 p)(Tyi ; )-Ixx) + (Ro^fQor) Jxi (2-97) 

In addition, for the auxiliary relations described by 
Equations 2-77, 2-78, 2-79, the corresponding perturbed 
equations are: 

p - 4> - TpsmOo - C^o cos ©o)9 (2-98) 

= 0 cos 6 0 - C V* sm 0 O sm <tl>)9 + cos 0 O cos d> 0 - 0„ sin <h> ) 4? +• 

+ 7 j> cos 0 O sir) d?o (2-99) 

r - t|>cos 0 O cos 4>o - C^P# cos0o sm d>» + 0 O cos d> 0 ) <p - 0 sm 4? 0 

■ C 4^0 sin ©o cos ) 0 (2-10 0) 

Although we have obtained linear equations, they are 
extremely complex and not suitable for control analysis. 

A further simplification is obtained if we consider the 
reference flight path to be straight, i.e. 

IP. - 0 O =0 (2-101) 

with no angular velocities 

- Qo * 1*0 = 0 (2-10 2) 

wings level flight, i.e. 

(p o =0 (2-103) 
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symmetric, i.e 





O 

ii 


(2-104) 




Vo * 0 


(2-105) 


With 


this assumption, the reference 


flight path is 


described 


by the equations: 






X„ = sin 0o 


(2-106) 




Y 0 = 0 


(2-107) 




lo - - 'WCjCOS 0 o 


(2-108) 




Lo = Mo - No — 0 


(2-109) 


while the 


perturbed equations become 






X = OM [ v + Wo + C| cos £1 0 J 


(2-110) 




Y = on [ v +l)<,r + VU 0 p - ^ cos Q> ] 


(2-111) 




2 = m [w - U> + cj sw 0 O 0 ] 


(2-112) 




L - Jx* p - Ixz r 


(2-113) 




M - 


(2-114) 




N = I Z2 f - 1x2 p 


(2-115) 


with auxilliary relations: 






p - (j> - ip sm 0» 


(2-116) 




r - 9 


(2-117) 




T r p cos 0 O 


(2-118) 



We find that the equations are separable into two groups. 
The first group associates the equations that describe the 
motion in the xz body plane, called longitudinal equations 
and the second group associates the equations that describe 
the motion in the xy body plane, called lateral directional 
equations. 87 



We find also, the force equations divided by mass and the 
moment equations divided by the corresponding moment of 
inertia. 

In the longitudinal equations the independent variables 

are: 



In the lateral-directional equations the independent 
variables are: 



with all these quantities being perturbations from the 
initial flight conditions. 

We arrive then to the following form of equations: 
Longitudinal equations 



v = x-axis linear velocity 
w = z-axis linear velocity 
q = pitching angular velocity 



v = y-axis linear velocity 
p = rolling angular velocity 
r = yawing angular velocity 



— X = v + Wo ^ t g cos 9 0 0 



( 2 - 119 ) 



±z = A -v. % + jsm&e 



(2-120) 






( 2 - 121 ) 



Lateral -directional equations 



— Y = v + H r - Wo p - g cos Go 4> 



( 2 - 122 ) 
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(2-123) 



Ixx K Ixx 



i I . Ixz • 

r^v *- = V Tw r 




(2-124) 



For further simplification, it is convenient to introduce 
here the wind axes system, which as we have seen is produced 
by a rotation about the y-body axis so that the x-wind axis 
points in the opposite direction of the relative wind. This 
axis is always tangent to the flight path and makes an angle 
of flight path with the horizontal, as shown in 
Figure 2-6. 



Figure 2-6 

Direction of the Body Axes with Respect to Wind Axes in 
Straight Level Perturbed Condition 




horizwrlal 



Since 




(2-125) 
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(2-126) 



in the steady situation 

00 = 

where 

cto =0 (2-127) 

In fact, in order to maintain straight level flight, 
has to be small. The appropriate relation for expressing the 
angle of attack <* in terms of the w-velocity if obtained by 
Equation 1-194 if one considers small c* and w , so that one 
can express: 

a - w/U„ (2-128) 

From Equation 2-127 then we deduce: 

<*o = Wo /lJo -0 =* Wo=0 (2-129) 

Notice that by this alignment, since the wind axes are 
fixed for fixed relative wind, we can measure the 
perturbation quantities about in the wind frame. 

It is further desirable to express the Lateral 

Directional equations in terms of the sideship angle £ 
instead of the v- velocity. The appropriate relation for 
expressing the angle {3 in terms of the v- velocity is 
obtained by Equation 1-195, if one considers small (3 and v , 
so that one can express 

§> = VU, (2-13 0) 

So introducing 

from Equations 2-126 and 2-127 

0O • j. (2-131) 
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from Equation 2-129 





W 0 = o 


(2-132) 


from Equation 2-130 








< 

H 

<F 


(2-133) 


we obtain 







Longitudinal Equations 

h X « i *■ gcosy.e 

s l - « - U.<j, +JS'njf.e (2-134) 




Lateral Directional Equations 

= v + Ur -geos#, 4> 



I XX 



, > Ivz 

L = f - — r 

1 Ixx 



(2-135) 





with auxilliary relations: 

p - ^ sm [o 

^ 0 (2-136) 

r = f cos ^ 

The only thing that remains is to express the left-hand 
side forces and moments in series of terms involving 
aerodynamic derivatives and small perturbation quantities, 
that in fact generate the forces and moments. 
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2.7 EXPANSION OF FORCES AND MOMENTS 



In Section 2.6, it was assumed that the motion of the 
vehicle consists of small deviations around a reference 
steady flight path. 

Small disturbance theory applied on the perturbed 
equations of motion assumes that the aerodynamic forces and 
moments are functions of the disturbance velocities (angular 
and linear), the control surface deflections and their 
derivatives . 

In general, for a function 

) (2-137) 

the Taylor series expansion about a reference value 

, .... ) (2-138) 

assuming that the function has derivatives of all orders 
defined, is written as 




(2-139) 



where the derivatives must be evaluated at the reference 
value denoted by the subscript (o) . 

Now, making the assumption that the deviations S3 , Sp,& 
Sx> , are so small that product terms are essentially zero, we 
obtain 




(2-140) 
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Using this method and provided that the above assumption 
is valid, we succeed to express the perturbation quantity SK 
as a linear combination of the perturbation parameters from 
which it depends upon. 

People in the field like to denote the partial 
derivatives as , Kp , Kv ", etc., instead of (3^/33 )° r 
(9K/3p ) a , ( OK /Qv ) a , etc. In addition, according to the 

convention met in Section 2.6 of dropping the $ , we obtain: 



where , K v , are called stability derivatives. 

Fortunately, not all the stability derivatives in an 
expansion of a force or a moment have a value. Some of them 
are zero or essentially zero. 

Expansion then of forces and moments, in terms of the 
significant stability derivatives yields: 



K - + K p jj K v v + 



(2-141) 



For longitudinal equations: 



TFT X - u + « + *s 5 



(2-142) 




(2-143) 



J 

M = + Mw w + w t 






(2-144) 



where it is clear that 




(2-145) 
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and 



denotes the control input due to control surfaces 



deflections. 

For the lateral -directional equations: 



— Y = Y v Y + Y; v + Y P p + Y r r + Ya 6 

on 


(2-146) 


^ L = 4 (S + Ljj £ + Lp P + Lr r v U S 


(2-147) 


~ N =• f Mp p + Nr r + M§5 

where it is clear that 


(2-148) 


y 4 *cY Y 4 *^Y 

v /m ’ P mi 9p 

Therefore the equations of motion become: 

Longitudinal Equations 


(2-149) 


'U + X v >u-X w w +|cos£,0 - Xgi 
- Z 0 u 4. Oj - Z w w - w - Vo ^ | sir? 0 - 25 § 

-M n w - Mi,w 4^ = M § 5 

Lateral Directional Equations 


(2-150) 


V - Yr V - Y* / - Yp p 4 u a r - Yr r - <j cos^ <j> - Y§ 5 

- 4 0 +p - 4P " Ix /f x , P - C r - US 

- N(? ^ - N(j |8 - p - Npp4-r-Nr' r - Ns <5 

with auxilliary relations of Equation 2-136. 


(2-151) 


It is more convenient to Laplace Transform the 


equations 



and arrange them in matrix form. Further augmenting the 
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equations by the auxilliary relations, using Equation 2-130 
for the lateral directional set and dividing the Y-equation 
by Vo we obtain the form 

Longitudinal Equations 



<J\ 

X 

c 


-Xw 


0 




s -sZw ~Zw 






-sM* -Mw 


S - M t 


0 


0 


i 



§ C0S k 




V 




X5 


<3 sm 




w 




Zs 


0 




% 




MS 


-s 




9 




0 



(2-152) 



Lateral -Directional Equations 



5- S Yv - Yv 


-Y* 

7 P 


4-Y* 


- s 4-v 


s-L p 




-s 


-*sMr 


s-M r 


0 


L 


0 


0 


0 





where 



i ccs fr 0 




$ 




Y/ 


0 0 




? 




is 


0 0 




r 


- 


NS 


-5 s(smp 




b 




0 


o -s(cos^o) 




Tp 




0 



£ (2-153 ) 



Y p * = Yp/Uo (2-15 4) 

Y* = Yr/l>o (2-155) 

Y* = Ys/Do (2-156) 

We will further simplify the equations of motion in the 
following chapters in a form that is more suitable for 
obtaining aerodynamic transfer functions, and computer 
solutions. 

Our interest first will be focused on examining the 
aerodynamic stability derivatives. They will be left in 
their dimensional form rather than non-dimensional. 
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CHAPTER 3 



NON-DIMENSIONAL STABILITY DERIVATIVES 
3.1 INTRODUCTION 

Stability is a very important concept in successful 
vehicle flight. We have encountered the terms "stability" 
and "equilibrium" but we did not emphasize particularly those 
concepts . 

In the last chapter we mathematically obtained equations 
of motion and linearized them about a referenced straight 
flight which was assumed to be an equilibrium condition. 

But what is equilibrium? Equilibrium is a state of rest 
or uniform (unaccelerated motion) in a straight line where 
all forces and moments acting on a body are balanced out. 

Stability is a property of an equilibrium state. We will 
consider two types of stability which are of great interest 
in flight vehicle dynamics, defined as follows: Static 
stability is defined by the direction and magnitude of the 
initial tendency to return to equilibrium condition after a 
small disturbance has occurred. It is distinguished by 
positive static stability if the body tends to return to its 
equilibrium state, negative static stability if the body 
diverges from its equilibrium state, and neutral static 
stability if the body remains in the disturbed position. 
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As an example, the pendulum considered in Section 1.3 has 
a positive static stability at its position but has a 

negative static stability at its ©=180" position, while its 
motion remains invariant if the pivot is moved in any other 
position and is said to have a neutral static stability. 

Dynamic stability determines the resulting motion in time 
of the body if initially disturbed. It is distinguished 
similarly in positive, negative, and neutral dynamic 
stability. 

Various conditions of dynamic behavior were considered in 
Chapter 1, in analyzing spring mass damper systems. 
Figure 3-1 summarizes all possible time responses of any 
dynamical system. 

Dynamic stability implies the existence of static 
stability, but the opposite is not necessarily true. 

Now the geometry of the vehicle comes to play a unique 
and important role in specifying any particular mode of 
dynamical motion. All surfaces exposed to the airstream 
contribute forces and moments about the center of gravity, 
which are expressed as we have seen in the last chapter as 
linear combinations of stability derivatives. 

The physical meaning of each of the stability derivatives 
is straight forward. Thus, X u is the rate of chance of the 
force generated in the x-direction due to a v-velocity 
perturbation, divided by the mass of the vehicle, while X y u 
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Z 

f 



-T5 



Ini4i«i disturbance 



a. brnniy damped 

_has_ 



c . 

r k 



s 



*a 

-o 



b. divergence 

*^iwf , 



Tvsiiivc da4ic- 1Wi*e dyimmic 

c. fleudral static rbbildy 
-Kme 



M«ja4ire slcrfic - Mecjcrbw dynamic 



Neutral static - Neutral dynamic 







IWitive static -Unit Ire dynamic 





%srH*e static- Neutral dynamic 

dirergent oscillation 
.'lime 



%siti»e static - Negatlre dynamic 



Figure 3-1 

Various Responses of a Dynamical System to a 
Small Disturbance 

a r b, c Non-Oscillatory Modes, d,e,f Oscillatory Modes 
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is the force generated in the x-direction due to a v-velocity 
perturbation divided by the mass of the vehicle. 

In this chapter we will give more insight in evaluating 
the stability derivatives in terms of the vehicle's geometry 
and airstream properties. 

3.2 LONGITUDINAL STABILITY DERIVATIVES 

Longitudinal stability derivatives are encountered in the 
longitudinal equations and, hance , they uniquely specify the 
X-force, the Z-force and the M-pitching moment. 

Dimensional analysis approach indicates that we can 
express the aerodynamic forces and moments in terms of their 
corresponding non-dimensional coefficients which are usually 
deduced from wind tunnel tests of models. 




Figure 3-2 

Forces and Moments Considered in Longitudinal Motion 
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In particular, we can express the longitudinal forces and 



moments as 



lift 




(3-1) 



drag 




(3-2) 



pitching moment M = Cm 2 P ^ ^ c 



(3-3) 



where C L ,Ci> ,C M are the corresponding lift, drag, and pitching 
moment coefficients, p the air density, V the velocity of the 
vehicle, S is a reference area usually taken the wing area, 
and c, a reference length, usually taken the length of the 
mean aerodynamic chord. 

The Longitudinal Stability Derivatives can be evaluated 
as follows: 

1. The Xu Derivative 



Is the rate of change in the X-force due to a 



v -velocity perturbation divided by the mass m 



. _ i 



(3-4) 



According to Figure 3-2, the X-force is expressed as: 



X - L sivi ot - j) cos +■ "C cos « T 



(3-5) 



where T denotes the thrust 



For at small, a T essentially zero, and making use of 



Equation 3-2 we obtain: 
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(3-6) 



x = -ip V^c„ »T 

Notice that because a small change only in v-velocity 

occurs 

V 2 - (Uo+v) 2 * V 2 + w 2 (3-7) 

and since products of perturbations are essentially zero 

V 2 = V 0 Z + 2Do v (3-8) 

where Uo is the reference flight forward velocity which for 
practical purposes is assumed equal to V. 

Then Equation 3-6 becomes 



X---ip(ttS2U»)SC t ,T (3 _ 9) 

Neglecting the effect of the v-velocity on the thrust 
T , by differentiating Equation 3-9 with respect to v, and 
dividing by the mass m, we obtain 



_ pUoS i ev 2 s % 

A ° cm ~ p 



2 m 'b v 



(3-10) 



Or since V=U„ 



(3-11) 



f Vo C Ko\ 

Xu " m * 2 ' 'bo ) 

Finally, if we define 

* 2 (3-12) 

we arrive at the following final form of the Xu derivative 



x„ 






(0> 'Cp„) 



(3-13) 
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The quantity Cj> u is essentially zero for low subsonic 
speed, but sometimes near the critical mach number 
(0.8<M<1.0) where a large increase in drag occurs, Cx>u 
obtains a large positive value. 

2. The Z 0 Derivative 

It is the rate of change in the Z-force, due to a 
u-velocity perturbation divided by the mass m. 



7 LUL 

Lv * (3-14) 

According to Figure 3-2, the Z-force is expressed as: 
Z = Wees© - L cosoc -'[Dsincx - Tsma T (3-15) 

which is simplified to: 

Z = W cos 0 - L (3-16) 

Neither W or 9 depend on ^-velocity . Since by 
similar approach as in the Xu derivative, we obtain: 



where 




/- Vo r <)Ci 

Lt > ~ T Hv 



(3-17) 

(3-18) 



The quantity Cc y is essentially zero for low subsonic 
speed, but sometimes near the critical mach number may reach 
a considerable value and its sign can be changed depending on 
the airframe geometry, mach number, dynamic pressure and 
aeroelastic properties. 

For subsonic straight level flight (&>=o) 
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pU/S/2 



(3-19) 



and Z v is expressed simply as 




(3-20) 



3. The M v Derivative 

It is the rate of change in the pitching moment M due 
to a ^-velocity perturbation, divided by the ly ^ moment of 
inertia. 

M --L351 

" ' I„ to (3-21) 

The Mo derivative is evaluated simply by 
differentiating Equation 3-3 with respect to o. 



of the u -velocity and its sign can change depending upon the 
geometry, aeroelastic properties, mach number and dynamic 
pressure. Thrust can also affect for propeller driven 
aircraft. 

4. The Xiaj Derivative 




(3-22) 



where 



/- Vo ^Cm 

" 7 (3-23) 

The quantity is very sensitive upon the variation 



(3-23) 



It is the rate of change in X-force, due to a 



w-velocity perturbation divided by the mass m 
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(3-24) 



y _ 1 3*. 
w ' w 'Jw 

A w-velocity perturbation causes a small change in 
the angle of attack <x . Thus, in this case we cannot neglect 
the angle of attach m which is approximately given by: 

oi Z 

V 0 ( 3 - 25 ) 

Therefore, the X-force is expressed as: 

X = L sine* -Dcos<* + T (3-26) 

or 



X = \ pV 2 £ (C L sm« - Q> cos o<) +T (3-27) 

Neglecting the effect of the w-velocity on the thrust 
T, by differentiating Equation 3-27 with respect to w and 
dividing by the mass m, we obtain 



v i 3x pU>s 

X w - — TT 7v“ - 



( [Cl +-0* sim +• Q> siruy - Cq v cosot ) 



m V 0 2n\ 

Now assuming <x small we express 



(3-28) 



Xw 



pUS 

2w 



(Cl-Cp„) 



(3-29) 



where 

_ OCp 

(3-30) 

The quantity C-p* represents the change in drag 
coefficient due to angle of attack changes. When the angle 
of attack increases, the total drag increases and this 
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quantity has normally a positive value. It is affected 
primarily by the wings, and a little by the tail and the 
fuselage. In aircraft dynamics it is not very important 
parameter mainly because the changes in angle of attack are 
small and, hence, its value is very small. 

5. The Z w Derivative 

It is the rate of change in Z-force due to a 
w-velocity perturbation divided by the mass m. 



z 



w 



1 2 L 

no 'Sw 



The Z-force simillarly is 

Z = W cos 0 - L cosm - D si no/ 



(3-31) 

(3-32) 



or 

Z = W CDS 9 - j pV 2 S Cq cosoc -Cj> smoi) ( 3 - 33 ) 

Neither W or 0 depend on w-velocity, since by 
differentiating Equation 3-33 with respect to w and dividing 
by the mass m we obtain: 



l 3Z pU>S r . 

.Csccs* ^ (3 _ 34) 



Now, assuming « small we express 



where 



Zw 



2m 





(3-35) 




(3-36) 
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The quantity C L</ represents the change in C L 
coefficient due to angle of attack changes. Graphically, it 
is well-known as "lift curve slope", as shown in Figure 3-3. 




Lift Characteristic Curves for Different AR's 

The changes in Q f are linear functions of oe , below 
the stall. Aspect ratio (AK), is defined simply as the ratio 
(span/chord) or in more general case as: 

A2 = (s P gn)2 = £ 

cvu>«3 area) g (3-37) 

This ratio, determines the slope of the lift curve. As shown 
in Figure 3-3, low A2 • s produce smaller lift coefficient G. 
for the same angle of attack ^ . 

The Q* derivative is always positive in the linear 
region. It is affected primarily by the wings, while the 
tail and fuselage contribute a little in aircrafts, but 
largely in missiles. Aeroelastic effects in addition may 
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cause bending or twisting of the wings that change the angle 
of attack and hence the value of . 

6. The Mw Derivative 

Is the rate of change in the pitching moment M, due 
to a w-velocity perturbation, divided by the Igy moment of 
inertia. 

A/I 1 

(3-38) 

A w-velocity perturbation is associated with an angle 
of attack which as we have seen is given by 

S w /Uo (3-39) 

When the angle of attack changes from its equilibrium 
position, the lift of the wing changes in the same sense, and 
that produces a positive or negative pitching moment 
depending on the center of gravity location. In addition, 
the lift of the tail changes and that contributes to the 
total pitching moment. These contributions together with the 
fuselage contribution make up this derivative. 

We desire a change in angle of attack to produce a 
pitching moment that tends to restore the original angle of 
attack. This requires a zero pitching moment in the 
undisturbed reference flight, a condition known as "trim 
condition" . 

The pitching moment, as we have seen, is expressed in 
terms of the non-dimentional pitching moment coefficient as 
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M = ipV ! S cC„ 



(3-40) 



Because a small change in w-velocity doesn't affect 
the V velocity, by differentiating Equation 3-40 with respect 
to w we obtain: 



_ QM Qqc _ i QCm 

<c)w ” 9* 9w ' ^ 'ScT ~ 2Vo ' C "97 (3-41) 

and finally by letting V-0 o we get 



0 S U> c 

Mvw = 5Iyy (3-42) 

The center of gravity location plays a very important 
role in derivative. In fact Cm# is proportional to the 

distance between the center of gravity and the aerodynamic 
center, the center where the lift force acts. 

If the center of gravity coincides with the 
aerodynamic center. Cm <* is zero; if it is ahead of the 
aerodynamic center, Cm# is negative and the airframe is said 
to be statically stable; and if it is aft of the aerodynamic 
center, Cm# is positive and the airframe is said to be 
statically unstable. The more ahead the center of gravity is 
located, the steeper the Cm# curve becomes. All the above 
cases are shown in Figure 3-4. 

It is clear that a stable situation like the center 
of gravity at B or A, an increase in angle of attack from 
to produces a negative pitching moment as Figure 3-4 shows, 
which tends to reduce the angle of attack. 
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Figure 3-4 

Cm. for Different Center of Gravity Locations 



Figure 3-4 also shows that for stable situations, for 
nearly small angles of attack, where the lift is nearly zero 
a mechanism must be provided to generate a positive . 
Here comes the effect of the tail which, being a negative 
angle of attack, generates a small lift force which though 
multiplied by a larger lever arm produces a balancing 
positive pitching moment. 

In other designs we find the tail being ahead of 
wings with positive angle of attack. This is a canard 
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configuration most usually found in missiles, also being the 
mechanism to generate a positive . 

In Figure 3-4 Cv,, curves are shown as linear 
functions of This is not true for large variations in c*, 

but can be nearly true for small variations around o/ ft . 

7. The Zm Derivative 

Is the rate of change in Z-force, due to a rate of 
change in w-velocity perturbation divided by the mass m 



The relation 






<*= w/Uo, 



1L 

' 3w (3-43) 

leads us to the expression 



- _L 51 

^ Vo 'So* (3-44) 

According to Equation 3-16 the Z-force is simplified 
to 

Z = W cos 0 - L (3-45) 

Neither W or 0 depend on rate of change in angle of 
attack cf , since differentiating Equation 3-45 with respect to 
of we obtain 

51 - -51 

0<x ” (3-46) 

where from Equation 3-1 



Therefore, 



51 = i p v 2 s 51 

o* * p ^ 






(3-47) 



7- J _ L __ i 3L -i i ^ L 

^ on c)^ on V 0 $<* mVo 'dot onU 0 2 P ^ (3-48) 
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To form a non-dimensional coefficient, we multiply 
and divide Equation 3-48 by c/2Uo so that letting V = Uo we 



get 








&c/2Vo) 



(3-49) 



or 




(3-50) 



where 




The quantity represents a change in lift 



coefficient due to a variation in the rate of change of the 
angle of attack « . 

If the angle of attack changes rapidly, the pressure 
distribution on the wing, tail and fuselage, does not adjust 
itself instantaneously to its equilibrium value. This, a 
change in lift and drag occurs to the sudden changes in angle 
of attack. 

Its sign is positive for low speeds but can be 
positive or negative on high speeds depending on aeroelastic 
effects . 

Also, a sudden change in angle of attack alters the 
downwash field as seen by the tail and a time delay is 
elapsed before the tail senses this alteration. 

All those effects are combined to give rise in the 
derivative. 
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8. The Mw Derivative 

Is the rate of change in the pitching moment M, due 
to a rate of change in w-velocity perturbation divided by the 
lyi moment of inertia. 

M . _ J_ W 

" lyy ' (3-52) 

This derivative W/Sw can be expressed as: 



Qw ” Vo' ry* (3-53) 

According to Equation 3-3 



M ~ ^ p V S c Civ^ 



(3-54) 



Since V is independent of the rate of change in angle 
of attack <*, by differentiating Equation 3-54 with respect to 
oc we obtain: 



21 = i pV 2 Sc 
Qi 2 r % 



Therefore, 



(3-55) 



M* = 2_.2! = _L_ . 21 = _L_ . 1 p V 7 s c 22? 

% U»Iyy *3* U>Ii/3 ^ (3-56) 

To form a non-dimensional coefficient, we multiply 
and divide Equation 3-56 by c/2Uo so that letting V = Uo, we 
get 

u . 1 I .1] ! C c 

M n — pUo OC — — • 

L^Iyy 2 r -?U> 0C*c/2V.) U 57 } 

or 



pSc 2 _ 

= — — C M « 

41^ 



(3-58) 



112 



where 



9Uc/2Uo) 



(3-59) 




For quantity Cm* represents a change in pitching 
moment coefficient due to a variation in the rate of change 
of the angle of attack «. 

The change in lift that is produced by sudden changes 
in angle of attack discussed already, multiplied by 
corresponding lever arms, generates changes in pitching 
moment and that is why Cm* arises. 

Same lag effects in the downwash field discussed, 
contribute to this derivative, as well as aeroelastic 
effects, that also disturb the air stream of the wing and 
tail, so that the pressure distribution over the aerodynamic 
surfaces does not adjust instantaneously. 

The sign of this derivative is negative and is 
relatively important in longitudinal dynamics. Large 
absolute values of this derivative are desired because it 
causes damping. 

9. The Zq Derivative 

Is the rate of change in Z-force, due to a q-pitching 
velocity perturbation divided by the mass m 



to: 



7 1 5L 

Z V n \ (3-60) 

According to Equation 3-16, the Z-force is simplified 



Z = W cos© -L (3-61) 
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Neither W or 9 depend on the pitching velocity (q) , 
since differentiating Equation 3-61 with respect to q, we 
obtain: 



n 

'H 



where from Equation 3-1 




Therefore: 



5L 

‘H 



• k P VIS - 



'h 



(3-62) 



(3-63) 



~7 - _ 1 _ ^ _ i I \jic '^~ t - 

~ mi ™ ~ ,w 2 9<^ (3-64) 

To form a non-dimensional coefficient we multiply and 

divide Equation 3-64 by c/2Uo so that letting V=Uo we get: 



or 



where 



l - .1 .1 p l) 0 z $ . 
% m 9 Y wo ^ on 



3Cc 



2U» yCcjc/lVo 



^ p S Uo C 

n = 




(3-65) 



(3-66) 



r %- 

L % ~ QC^c/2a) (3-67) 

The quantity represents a change in lift 

coefficient due to a variation in the pitching velocity. 

If the airframe pitches about its center of gravity, 
the angle of attack changes and a change in left of wing, 
tail and fuselage is developed that gives rise to this 
derivative. 
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If has a small positive value in low speed flight but 
may have large values positive or negative in high speeds due 



to aeroelastic effects. 

Curved longitudinal flight paths give rise also to 
centrifugal forces on all components of the airframe that due 
to aeroelastic effects cause a bending of aerodynamic 
surfaces and, hence, a change in lift. 

In the past without computer solutions available, 
this derivative was considered zero to simplify the analysis. 

10. The Mq Derivative 

Is the rate of change in the pitching moment M, due 
to a pitching velocity perturbation, divided by the Iyy 
moment of inertia. 



Since V is independent of the pitching velocity q, by 
differentiating Equation 3-69 with respect to q, we obtain: 




(3-68) 



According to Equation 3-3, 



M i p V 2 5cC M 



(3-69) 




(3-70) 



Therefore 




(3-71) 
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To form a non-dimensional coefficient, we multiply 
and divide Equation 3-71 by c/2Uo, so that letting V=Uo we 
get: 



or 



where 



m,* 4- pa 2 sc — 



M 






pV 0 $c 2 

4 % 




(3-72) 



(3-73) 



C — 

M r ^C^c/2a) (3-74) 

The quantity Cm^ represents a change in pitching 
moment coefficient due to a variation in the pitching 
velocity q. 

The change in lift that was produced, due to pitching 
velocity perturbations discussed already, multiplied by 
corresponding lever arms, generates changes in pitching 
moment and that is why C rises. 

The tail is the main contributor to this derivative 
which is also referred to as "pitching damping derivative". 
Aeroelastic effects and centrigugal forces generated due to 
curved longitudinal flight paths also contribute to this 
derivative . 

In longitudinal dynamics is very important. It 
has a negative value, and the larger it is, the more the 
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damping is increased; a condition that is highly desired. In 
high speed flights may have positive or negative value 
depending on aeroelastic effects. It is also noted that 
tailless airplanes have a poor damping. 

3.3 LATERAL - DIRECTIONAL STABILITY DERIVATIVES 

Lateral -Directional stability derivatives are encountered 
in the lateral directional equations and, hence, they 
uniquely specify the Y-force, the L-rolling moment and the 
N- yawing moment. 

Dimensional analysis approach indicates that we can 
express the lateral-directional aerodynamic forces and 
moments in terms of their corresponding non-dimensional 
coefficients as: 



Y-force 


Y 




(3-75) 


L-rolling moment 


L 


= c g { pv 2 s y 


(3-76) 


N- yawing moment 


N 


= C„ i pV*Sb 


(3-77) 


Cy , Q , Q are 


the 


corresponding Y-force, 


rolling 



moment and yawing moment coefficients, and b is a reference 
length, usually taken the wing span. 

The Lateral-Directional Stability Derivatives can be 
evaluated as follows: 

1. The Yv Derivative 

Is the rate of change in Y-force due to a v-velocity 
perturbation divided by the mass m: 
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(3-78) 



Y * ?L 

v = m 9v 

A v-velocity perturbation causes a small change in 
the angle of sideslip £ which is expressed according to 
Equation 2-133 , approximately as: 

S v/Vo (3-79) 

so that 

Y _ •* r Si’ 

“ JwuT‘ QjT (3-80) 

Assuming that the V-velocity is independent of the 
angle of sideslip^, differentiating Equation 3-75 with 
respect to £ , we obtain: 

2L , i p v 2 s 

•)(! 2 r 9|S (3-81) 

Therefore, letting V=Uo we get: 

y P v °$ r 

lv ' 2m Hs (3-82) 

where 

(3-83) 

The quantity represents the change in the Gf side 

force coefficient due to a change in the angle of sideslip f. 

The main contributor of this derivative is the 
vertical tail. Secondarily the fuselage, and also the wings 
contribute a little. 
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Figure 3-5 shows the forces that the airframe 
experiences due to a v-velocity perturbation. First an angle 
of sideslip builds up and all aerodynamic surfaces are 
sensing a change in the approaching free stream. 







Figure 3-5 

Forces Experienced Due to a v-Velocity Perturbation 

As a consequence, a force F^ due to the fuselage and 
a force F v t due to the vertical tail appears. The wings also 
experience a different lift force. For the configuration 
shown in Figure 3-5, the right wing experiences more lift 
than the left wing because it senses a greater V-normal 
component of the approaching free stream. Those lift forces 
may contribute to a small change in the Y-force. 

All these contributions make a rise to derivative 
which is very important in lateral-directional dynamics. It 
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introduces a damping that tends to eliminate the angle $ , 
because the resultant Y-force appears acting on a point aft 
of CG location and that generates the appropriate moment to 
do this. 

Its value is negative and it might be reasonable to 
desire a large value , for the £ angle to die out quickly. 
Although high maneuverability airplanes must have a small 
negative value for this derivative, if they were to perform 
easy roll maneuvers and banked turns. 

2. The Ly Derivative 

Is the rate of change in L-rolling moment due to a 
v-velocity perturbation divided by the Ixx moment of inertia. 



Lv s lx*' 9 y (3-84) 

It is also clearly expressed as: 

I . J_.2L 

v 'tU«* 3(2 (3-85) 

Differentiating Equation 3-76, with respect to /£ , we 



obtain 



.2L _L n V 2 ^ L 2£l 

^ : 2 P V b 3(5 (3-86) 

Therefore, letting V=Uo we get: 

, past 

iv = % (3-87) 
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where 



C ^ 

h " ^ (3-88) 

The quantity represents the change in the Cg 

rolling moment coefficient due to a change in the angle of 
sideslip $ . 

It is usually referred to as the effective Dihedral 
Derivative. 

Dihedral is often used as a mechanism to improve 
lateral-directional stability. Figure 3-6(a) shows a head-on 
view of an airplane whose wings are tilted up at some 
dihedral angle with respect to the horizontal. 

We have already discussed that due to a v-velocity 
perturbation, according to the configuration shown in 
Figure 3-5, where a positive sideslip angle is shown, right 
wing generates more lift than the left. This is true for the 
swept wings. If the wings are not swept, the mechanism to 
produce the same effect is the dihedral. Usually even swept 
wings are given a small dihedral to improve 
lateral-directional stability. 

To understand the dihedral effect, consider the 
aircraft shown in Figure 3-7 with non-swept wings but with 
some dihedral angle. The aircraft is distributed by an angle 
of sideslip ^ . 
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Figure 3-6 

Dihedral Effect on Lateral-Directional Stability 




Dihedral Effect on Non-Swept Wings 
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The perpendicular components of the free stream 
velocity are equal but the angle of attack that the right 
wing projects to this perpendicular component is greater than 
the angle of attack that the left wing projects, as an effect 
of the dihedral angle. 

As a consequence, the right wing generates more lift 
than the left wing as shown in Figure 3-6(c). Those lift 
forces contribute to a change in the rolling moment 
coefficient. 

Figure 3—6 ( d ) and (c) show the effect of the 
high-wing placement which has the same effect as the dihedral 
angle. It is clearly shown that in the low wing placement 
configuration, the fuselage does not permit the vertical 
component of the free stream to act upon the right wing to 
produce a stabilizing moment, whereas in the high wing 
placement configuration the fuselage effect is nearly zero 
and that makes a stabilizing rolling moment to show up. 

In addition, the force acting on the vertical tail 
shown in Figure 3-5, multiplied by the corresponding vertical 
arm relative to the x-longitudinal axis, contributes also to 
the same effect. 

typ is nearly always negative in sign and the 
corresponding rolling moment L, is negative for a positive 
sideslip angle ft , It is a very important derivative in 
lateral directional dynamics because it causes a damping and 
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improves the stability. Although high negative values of 
this derivative will tend to decrease the maneuverability of 
the airplane in roll movements. In missiles a low negative 
value of this derivative is desired. 

3. The Nv Derivative 

Is the rate of change in N-yawing moment due to a 
v-velocity perturbation divided by the Izz moment of inertia. 



N J_ 
v 3 Jzz 9v 

It is also clearly expressed as: 



(3-89) 



> j)W 

VoJzZ 



(3-90) 



Differentiating Equation 3-77 r with respect to we 



obtain : 







Therefore, letting V=Uo we get: 



(3-91) 



P ast, 

- IST Ct, |3 (3-92) 

where 

^ /90 (3-93) 

The quantity Cn» represents the change in the Cn 
yawing moment coefficient due to a change in the angle of 
sideslip 
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It is usually referred to as the Weathercock 
Derivative. The main contributor of this derivative is the 
vertical tail and secondarily, the fuselage. If the airframe 
experiences a v-velocity perturbation and an angle of 
sideslip appears, the free relative stream strikes the 
fuselage and the vertical tail in such a manner that the 
force Ff , shown in Figure 3-5, produces a destabilizing 
yawing moment and the force *>i a stabilizing moment. Wings 
contribute a little in the positive (stabilizing) sense. The 
total value of this derivative is positive, signifying static 
directional stability. 

The Cvp derivative, is very important in 
lateral-directional dynamics. The vertical tail must be 
properly designed to provide a sufficient weathercock effect. 
At large, sideslip angles a stall may occur due to flow 
separation and a catastrophic sideslip divergence may result. 
A low aspect ratio prevents stalling but decreases stability. 

Usually a large positive value of the ^ derivative 
is desired, by extending the vertical tail. A large value 
helps the pilot in good coordinated turns and prevents 
sideslip. 

4. The Yy Derivative 

Is the rate of change in Y-force, due to a rate of 
change in v-velocity perturbation, divided by the mass m. 




21 

"3v 



(3-94) 
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It is also clearly expressed as: 



y. J_ j>Y_ 

' ' l>„w 9(S (3-95) 

Differentiating Equation 3-75 with respect to ft we 



obtain: 






oT = 7 P V ' S 



*3 



2 r ' ~ % (3-96) 

To form a non-dimensional coefficient we multiply and 
divide Equation 3-96 by b/2Uo. 



^Y 






il, i p v*si- . 

2 f 21X 9C|Jb/2&) 
Therefore, letting V=Uo we get: 



(3-97) 



Y. = pSb Ci 
* 4nn ^ 



(3-98) 



where 



r - 

% ' QC(jb/20,) (3-99) 

The quantity represents the change in the 

coefficient with variation in rate of change of the sideslip 
angle (3. 

If the sideslip angle changes rapidly, the pressure 
distribution on the vertical tail and fuselage does not 
adjust itself instantaneously to its equilibrium value. 

This effect, as well as aeroelastic effects, give 
rise in derivative but its value is very small and is 

usually neglected. 
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5. The Lv Derivative 



Is the rate of change in L-rolling moment, due to a 
rate of change in v-velocity perturbation divided by the Ixx 
moment of inertia. 



1 . _ 1 9L 

Ixx 9v 

It is also clearly expressed as: 


(3-100) 


L . _ i 3L 

v “ a ixx 


(3-101) 



Differentiating Equation 3-76, with respect to /3, we 

obtain: 



* .ipVSbi- 
% 2 p % 

To form a non-dimensional coefficient 


(3-102) 


% t y2 C | b 

Of ’ 2 P ^ 2Vo 0(fb/2V.) 

Therefore, letting V=Uo 


(3-103) 


Lv ' 4lx* 

where 


(3-104) 


r % 

% 0C(ib/2Vo ) 


(3-105) 



The quantity represents the change in the Q 
rolling moment coefficient with variation in rate of change 
of the sideslip angle 
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The change in the side forces on the vertical tail 
and the fuselage produced due to sudden changes in the 
sideslip angle /3 , discussed already, multiplied by 
corresponding lever arms provide a total change in the Q 
rolling moment coefficient. 

This derivative has a very small value because the 
lever arms are very small and it is usually neglected. 

6. The Ny Derivative 

Is the rate of change in the N-yawing moment, due to 
a rate of change in v-velocity perturbation divided by the 
Izz moment of inertia. 



Nv = 



L 

Izl Av 



(3-106) 



It is also clearly expressed as: 



N v - 



1 



SN 



V 0 1i2 (3—107) 

Differentiating Equation 3-77, with respect to /3, we 



obtain : 



ON < ./i c , OCn 

- 7 - = — p V S b 777 - 

To form a non-dimensional coefficient: 



(3-108) 



M = i p V 2 Sb — 

% 2 r OIL 



QCr, 



2 r ' ' QVc QC/ib/2a) 

Therefore, letting V=Uo we get: 



m- P sbi r 

N ' • HIT C 



(3-109) 



(3-110) 
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where 



r _ 

* 9C/2b /2tt) (3-111) 

The quantity represents the change in the Gi 

yawing moment coefficient with variation in rate of change of 
the sideslip angle . 

The change in the side forces on the vertical tail 
and the fuselage produced due to sudden changes in the 
sideslip angle @ , discussed already, multiplied by 
corresponding lever arms provide a total change in the Gt 
yawing moment coefficient. 

The lever arm of the vertical tail this time is long 
and hence, this derivative may be considerable. Positive 
values of this derivative increase the damping in the yawing 
motion and are desirable in general. 

7. The Y ? Derivative 

Is the rate of change in the Y-force, due to a 
p-rolling velocity perturbation divided by the mass m. 



v _ L 51 

T P ' m 3p (3-112) 

Differentiating Equation 3-75, with respect to p, we 



obtain : 



C )Y 





(3-113) 
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To form a non-dimensional coefficient 



Hi _ i d V 2 S -L - 

9p ~ 2 r 2 Uo 9(pb/2U» ) (3-114) 

Therefore, letting V=Uo 

Y pU>St 

? = A'M c yp (3-ii5) 

where 

c - 

^ ' 9Cf>b/2U.) (3-116) 

The quantity Cy p represents the change in the Cy 
Y-force coefficient with variation in the rolling velocity p. 

When an aircraft has a roll velocity p, a force is 
acting on the vertical tail that usually opposes the roll 
movement. This consists the main contribution to this 
derivative and is small enough to be neglected. 

The force acting on the vertical fin may be affected 
by the sidewash from the wing that gives a slightly different 
direction of the airstream. 

8. The Lp Derivative 

Is the rate of change in the L-rolling moment, due to 
a rolling velocity p perturbation divided by the Ixx moment 
of inertia. 



I J_ 
p ‘ ixx 9 p 



( 3 - 117 ) 
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Differentiating Equation 3-76 , with respect to p, we 



obtain: 



= i pV'Sb ^ 

On 2 > 



QL 

Op 2 r ~ ~ Op 
To form a non-dimensional coefficient. 



(3-118) 



5 L 

Of 




b OQ 
' QCpb/2(I) 



Therefore, letting V=Uo 



p * m** 




(3-119) 



(3-120) 



where 



C» - — ^ — 

p 0Cpb/2U,) (3-121) 

The quantity Q p represents the change in Q rolling 
moment coefficient with variation in the rolling velocity p. 

This derivative is also referred to as "roll damping 
derivative", because it has usually a negative or damping 
effect coming mainly from the wings and secondarily, from the 
vertical and horizontal tails. 

When the aircraft rolls, the force acting on the 
vertical tail, shown in Figure 3-8, multiplied by the 
corresponding lever arm, produces an opposing rolling moment. 
In addition, on the down-going wing, the effective angle of 
attack is increased and on the up-going wing, decreased. 
This makes the lift on the down-going wing greater and thus 
an opposing rolling moment is produced. Similar effect, on a 
smaller scale, is produced by the horizontal tails. 
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T 




Figure 3-8 

Forces Due to Rolling Velocity (p) 

This is true as long as the down-going wing is not 
flying near the stall due to a large angle of attack 
experienced in some cases. In such cases a positive effect 
will result in autorotation. 

The derivative Qp is quite important in 
lateral-directional dynamics. The value of this derivative 
directly affects the ailerons design. A small negative value 
of this derivative will result in a better aileron input 
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response while a large negative value will make the airframe 
heavily damped in roll maneuvers. 

9. The Np Derivative 

Is the rate of change in the N-yawing moment, due to 
a rolling velocity -p perturbation, divided by the Izz moment 
of inertia. 



N 



P 



J_ 3N 

I 2 Z 9p 



(3-122) 



Differentiating Equation 3-77 , with respect to p, we 



obtain: 



2L , IpV'Sb 

Op 2 r 9 f (3-123) 

To form a non-dimensional coefficient, 

M _ j_ y2g i OGt 

= P 20, ' 30pb/2UO (3-124) 

Therefore, letting V=Uo 

N - ftt,sb * r 

S’ (3-125) 

where 

__ QCu 

~ OCpb 72v7) (3-126) 

The quantity Cn p represents the change in C« yawing 
moment coefficient with variation in the foiling velocity p. 

This derivative usually has a negative or positive 
effect coming from the wings, the vertical tail and 
secondarily, from the horizontal tail. 
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When the aircraft rolls, the force acting on the 
vertical tail, shown in Figure 3-8, multiplied by the 
corresponding lever arm produces a yawing moment towards the 
same side of the roll motion, i.e., roll to starboard causes 
yaw to starboard. This contributes a positive effect and 
depands upon the vertical tail geometry and the sidewash from 
the wings. 

In addition, the contribution of the wings can be 
considered to be made of two parts, one opposite the other, 
the total effect being negative. When the airframe rolls, on 
the down-going wing, the effective angle of attack increases 
and on the up-going wing decreases, as shown in Figure 3-8. 
As a consequence, the lift and drag are greater on the 
down-going wing than on the up-going wing. The projections 
of the drag forces on the xcy plane usually produce a yawing 
moment towards the same side of the roll motion. This 
produces a positive effect while the projections of the lift 
forces produce a yawing moment towards the opposite side of 
the roll motion, and this produces a negative or damping 
effect. 

The overall value of this derivative is usually 
negative, i.e. roll to starboard produces a yaw to port. At 
higher speeds the numerical value decreases and may become 
positive . 
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This derivative is fairly important in 
lateral-directional dynamics, as far as roll maneuvers and 
their effect to yaw. 

10. The Yr Derivative 

Is the rate of change in the Y-side force, due to a 
yawing velocity perturbation r, divided by the mass m. 



Y - L 

r = <m ' Qr (3-127) 

Differentiating Equation 3-75, with respect to r, we 

obtain: 



2I.-i P v 2 s 

■Jr ! r » 



To form a non-dimensional coefficient. 



i|Y _ J_ 0 V 2 S A ^ 

9r 2 r 2Vo ' 30b/2a) 

Therefore, letting V=Uo 



where 



v pUoSb 

r = ~ 4*~” 



Cy r 



(3-128) 



(3-129) 



(3-13 0) 



r ^ 

^ r ' QCrk/20.) (3-131) 

The quantity represents the change in Qj Y-force 

coefficient with variation in the yawing velocity r. 

When the aircraft has a yaw velocity r, as Figure 3-9 
shows, a force is acting on the fuselage that usually has a 
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negative or damping effect on the yaw motion and a force is 
acting on the vertical tail that usually has the same effect. 




Figure 3-9 

Forces Due to Yaw Velocity (r) 

Both contributions are small enough to be neglected, 
therefore, this derivative is of little importance in 
lateral-directional dynamics. 

11. The t* r Derivative 

Is the rate of change in the L-rolling moment, due to 
a yawing velocity r perturbation, divided by the Ixx moment 
of inertia. 

i _L A 

L|r ' I** ' (3-132) 
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Differentiating Equation 3-76, with respect to r, we 



obtain : 



— = i pV 2 Sb 

9r 2 r Qr 

To form a non-dimensional coefficient, 



(3-133) 



_L 

Or ' 2 P 2^o ' OCrb/ZUO 

Therefore, letting V=Uo 



lr 



p Vq S b 2 

4lxx _ e r 



where 



(3-134) 



(3-135) 



_ 2£i 

" e>C rfc> / 2L>o ) (3-136) 

The quantity Qr represents the change in Q rolling 
moment coefficient with a variation in the yawing velocity r. 

When the aircraft is yawing about its vertical axis, 
the outer wing moves faster than the inner. As a 
consequence, the lift on the outer wing is greater than that 
of the inner and, hence, a rolling moment is produced towards 
the same side of the yaw motion. This contributes a positive 
effect, i.e. yaw to starboard produces roll to starboard. 

In addition to this major wing contribution, the 
force acting on the vertical tail multiplied by the 
corresponding lever arm, produces the same rolling effect. 
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The total value of this derivative is usually 
positive, but it is not so important in lateral -directional 
dynamics. 

12. The Nr Derivative 

Is the rate of change in the N-yawing moment, due to 
a yawing velocity r perturbation, divided by the Izz moment 
of inertia. 



Nr -- 



ON 



Izz ' 3r (3-137) 

Differentiating Equation 3-77, with respect to r, we 



obtain : 



; i pV'Sb 

Or 2 ' 9r 

To form a non-dimensional coefficient. 



™ * i p v J sb i. . %K — 

Or 2 ' 211 c )Crh>/2Vo ) 

Therefore, letting V=Uo 

M . P».Sb ! 

Nr ’ "1ST Civ 



where 



Cn r - 



OCn 



(3-138) 



(3-139) 



(3-140) 



OOb/2UO (3-141) 

The quantity Gr\ r represents the change in yawing 
moment coefficient, with a variation in the yawing 
velocity r. 
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This derivative is known as the "yaw damping 
derivative", because it usually has a negative effect. 

The main contributor is the vertical tail. When the 
aircraft is yawing, the force acting on the vertical tail 
multiplied by the corresponding long lever arm, produces a 
considerable yawing moment that opposes the initial yawing 
disturbance . 

In addition, the drag on the outer wing is greater 
than that of the inner wing. That produces a yawing moment 
that opposes the initial yawing disturbance. The force on 
the fuselage also contributes an opposing effect. 

This derivative is very important in 
lateral-directional dynamics. Its total effect is negative 
and large values are desired for more effective damping. 

3.4 CONTROL SURFACE DEFLECTIONS 

We have discussed so far all the important stability 
derivatives that rise due to the airframe geometry. They 
uniquely specify the forces and moments about the center of 
gravity provided that they are multiplied by the correct 
factor that depends upon the airstream and airframe 
characteristics. 

A stability analysis is possible now for determining the 
longitudinal and lateral directional modes of motion, 
corresponding to an initial disturbance, i.e. homogeneous 
solution to the differential equations. 
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In addition, some means of control in longitudinal and 
lateral directional motion are available to the human or 
automatic pilot. They are shown in Figure 3-10. 






Figure 3-10 

Basic Longitudinal and Lateral Control Systems 
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Elevator, controls the longitudinal motion of the 
airframe providing nose up or nose down moments as a response 
to corresponding control surface deflections. 

Aileron and rudder, control the lateral directional 
motion providing rolling and yawing moments correspondingly 
as a response to the corresponding control surfaces 
deflections. 

The symbols used for the corresponding angles of rotation 
and control surface deflections are shown below: 



axis 


Ox 


0 y 


Oz 


anglle of rotation 


<t> 

C bank) 


0 

C pitch) 


ip 

( yaw) 


control surface defP. 

- 


5 

C ai leron) 


n 

(elevator) 


5 

C rudder) 



In the longitudinal and lateral directional equations of 
motion. Equations 2-152 and 2-153, respectively, we have 
denoted by the symbol $ the control input due to control 
surfaces deflections. 

The right-hand side of the longitudinal Equation 2-152 is 
expanded as follows: 

*4 - X* 

Z S = Z« (3-142) 

= Mn + Mti 



141 



and the right-hand side of the lateral directional 
Equation 2-153 as: 



YS = Y 5 

lh = L§ + Lg + L S (3-143) 

Nj = Ng 4. Ng + Ng 

In addition, to the right-hand side terms that take care 
of the main control surfaces, there may exist terms rising 
from the auxiliary surfaces such as flaps, stabilizers, dive 
brakes, etc. 

A general approach suggests to examine the subscripted 
derivatives instead of the particular derivatives dedicated 
for a specific control surface deflection, i.e. the total 
effect by all the possible combinations of the control 
surface deflections. The approach is similar to that 
followed in Sections 3.2 and 3.3. 

1. The Derivative 

Is the rate of change in the x-force due to a 
control surface perturbation divided by the mass m. 



= 






m 'BS (3-144) 

By § we denote a general control surface deflection. 
For longitudinal control it might be elevator, stabilizer, 
flap, slat, dive brakes, etc. 

According to Equation 3-27, the X-force is simplified 

to: 

X = - J P v S Cp +T (3-145) 
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We can assume that neither V or T depends on the 
control deflection. Thus differentiating Equation 3-145, 
with respect to , we obtain: 



9X 

96 



4* pv 2 s 



96 



If we let V=Uo and 

= / 9<5 

we arrive at the expression: 



(3-146) 

(3-147) 



Xl - - c 

S 2« Cp5 (3-148) 

The quantity represents a change in the drag 

coefficient due to a variation in a longitudinal control 
surface deflection. 

A control deflection, positive or negative, increases 
the drag due to the greater projected frontal area, but this 
sometimes is desirable because it provides damping in an 
oscillatory perturbed motion. 

2. The lb Derivative 

Is the rate of change in the Z-force due to a S 
control surface perturbation divided by the mass m. 



to: 



7 _ 9Z_ 

^6 " m ' 

According to Equation 3-16 , 



(3-149) 

the Z-force is simplified 



Z r W cos 0 - L 



(3-150) 



143 



or 

Z = v\y cos® - j p V 2 s c L 



(3-151) 



Thus di 
obtain: 



Neither W or 9 depends on the 
f fe r entiating Equation 3-151, 



- i_ oV 2 S — 

28 2 P 3S 



^ control deflection, 
with respect to 5 we 

(3-152) 



If we let V=Uo and 

Qj = e)Ct/05 (3-153) 

we arrive at the expression 



7 -L^£r 

3 ' S'* LL5 (3-154) 

The quantity Qj represents a change in the lift 
coefficient due to a variation in a longitudinal control 
surface deflection. 

This derivative is a measure of the effectiveness of 
the elevator mainly, in producing a change in the lift force 
acting on the tailplane. A positive elevator deflection 
produces an increase of lift, i.e. negative Z-force. 

It has a very small value on aircrafts with tail but 
a large value on tailless aircrafts. 

3. The M § Derivative 

Is the rate of change in the M-pitching moment due to 
a control surface perturbation divided by the Iyy moment of 
inertia 



M b 



[ 



(3-155) 
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According to Equation 3-3, the pitching moment is 
expressed as: 

M = j p V 2 Sc C M (3-156) 

Differentiating Equation 3-156, with respect to $ we 

obtain: 

'cM w2<- 

VS ~ 2 P C 9J (3-157) 

If we let V=Uo and 

Cm 5 ^/3<5 (3-158) 

we arrive at the expression 

pU„ 2 $C 

23yy CM5 (3-159) 

The quantity represents a change in the pitching 

moment coefficient due to a variation in a longitudinal 
control surface deflection. 

It is commonly referred to as the elevator power 
derivative because it is a measure of the elevator control 
effectiveness. 

The change in Z-force produced by the elevator 
control surface deflection, multiplied by a rather long arm, 
provides a pitching moment in the negative sense, i.e. a 
positive elevator deflection results in a nose-down 
(negative) pitching moment. 

That is the primary function of the elevator and is 
considered to be the most important of all the control 
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surface functions. It adjusts the angle of attack and it 
trims the vehicle in its 1 -ongitudinal motion, by providing 
opposite pitching moments to accommodate disturbing moments 
and center of gravity relocations. 

For a design point of view, the larger the center of 
gravity range is, the greater the elevator pitching moments 
must be provided, i.e. larger G 45 values. 

4. The Ys Derivative 

Is the change in the Y-side force due to a 2 control 
surface perturbation divided by the mass, m. 



Y 1 
15 ' m ’as 

Differentiating Equation 3-75, 



(3-160) 

with respect to S we 



obtain: 



51 




QCy 

2)5 



If we let V=0o and 

Ctj 5 c 2)G, /9§ 

we arrive at the expression: 



(3-161) 

(3-162) 



Ys _pE£- cm . ,, 

2 om (3-163) 

The quantity represents the rate of change in the 

Y-side force coefficient due to a variation in a lateral 

control surface deflection. 
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This derivative rises mainly from the rudder. A 
positive rudder deflection is defined to produce a positive 
side force. Therefore, the value of this derivative is 
positive but small enough to be neglected. 

Rudder deflections also cause an increase of the 
total drag of the vehicle but that has a negligible effect on 
the longitudinal motion, and is neglected. 

5. The U Derivative 

Is the rate of change in the L-rolling moment due to 
a control surface perturbation, divided by the Ixx moment 
of inertia. * 

J_ 3l_ 

Ixx (3-164) 

Differentiating Equation 3-76 , with respect to 5, we 

obtain : 




5k 1 

* 2 

If we let V=Uo and 




c e 5 - 

we arrive at the expression: 



(3-165) 

(3-166) 



pVoSb . 

lh r 2 Ixx Ce<J (3-167) 

The quantity Qj represents the change in the Q 
rolling moment coefficient due to a variation in a lateral 
control surface deflection. 
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The side force generated by a rudder deflection, 
produces a rolling moment in the positive sense, i.e. a 
positive rudder deflection generates a positive rolling 
moment. In airplanes, with large vertical tails, this effect 
is negligible and is often neglected. 

The main contribution to this derivative is from the 
ailerons, and it is a measure of the aileron effectiveness. 
A positive aileron deflection is defined as the one to 
produce a positive rolling moment, i.e. left aileron down and 
right aileron up, therefore, the value of this derivative is 
positive. 

In lateral-directional dynamics, this derivative is 
probably the most important because it establishes the 
maximum available rate of roll of an airplane which is a very 
important consideration in combat environment. 

It also plays an important role in take-offs, and 
landings, where adequate and rapid lateral control is highly 
desirable. 

Rapid rate of change of aileron deflections gives 
also a rise to this derivative due to lag effects, i.e. not 
immediate adjustment of the airstream. 

6. The M6 Derivative 

Is the rate of change in the N-yawing moment due to a 
£ control surface perturbation, divided by the Izz moment of 
inertia. 
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Differentiating Equation 3-77, with respect to S , W e 



obtain : 

4 \ti c < ^' T ’ 

ft" " 2 P bC W (3-169) 

If we let V=Uo 

Crj 5 =3Cn/ as (3-170) 

we arrive at the expression: 

2 Izz (3-171) 

The quantity Cflj represents the change in the C n 

yawing moment coefficient due to a variation in a lateral 
control surface deflection. 

The side force generated by a rudder deflection, 
produces a yawing moment in the negative sense, i.e. a 
positive or to the left rudder deflection, produces a 
negative yawing moment. This derivative is a measure of the 
rudder effectiveness. 

Rapid rate of change of rudder deflections gives also 
a rise to this derivative due to lag effects. 

This derivative plays an important role in crosswind 
takeoffs, and landings, and counteracts adverse yaw in 
rolling maneuvers. 




Finally, aileron deflections give also a rise to this 
derivative since each aileron introduces a drag force of 
different magnitude. This effect may be positive or negative 
depanding on the rigging of the aileron and the angle of 
attack of the airframe. 

3.5 DETERMINATION OF THE STABILITY DERIVATIVES AND OTHER 
CONSIDERATIONS 

All the stability derivatives discussed so far rise from 
the flight vehicle's geometry, being highly dependent upon 
the design layout. We have referred mostly in this chapter 
to the conventional airplane with horizontal tail and we have 
seen that the derivatives are made up of a number of separate 
contributions. 

We had strictly remained in the dimensional form of the 
stability derivatives and we simply expressed them in terms 
of airstream characteristics, airframe geometry and 
non-dimensional derivatives. 

Dimensional form has the advantage of obtaining real time 
solutions rather than working with normalized time mass and 
inertia components that must also be non-dimensional in 
considering that form. 

The advantage of non-dimensional derivatives, i.e. 
independent of size and airstream characteristics is 
maintained because dimensional derivatives are expressed with 
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the corresponding non-dimensionals, some of which though do 
depend upon the airstream and airframe characteristics. 

Thus, some of the dimensional expressions and some of the 
non-dimensionals are constants, but they can be expressed in 
polynomial form. 

As an example, for the METEOR 7: 

Q v = -0.052 - 0.009 C L 
Cg r = 0.0295 + 0.159 C L 
Grip — —0.042 C L 
Gn ? = -0.182 - 0.052 C* 

For others, plots may exist indicating their variation 
with respect to C L or mach number. 

During the development of an aircraft, it is important to 
know the numerical values of the derivatives under various 
conditions of flight. In the first stage, theoretical 
estimates can be made supported by wind tunnel tests and 
model experiments. 

By then, simulation in digital computers can be made to 
obtain response characteristics to initial disturbances or 
control deflections and determine whether the values of the 
derivatives are satisfactory. 

When the aircraft has been built, flight tests are 
carried out to confirm the theoretical with the actual 
responses. 
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There also exist flight test techniques for determining 
the values of the derivatives. 

Tables 3.1 and 3.2 summarize the longitudinal and lateral 
directional derivatives including definition, unit, 
expression and auxilliary relations. 
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TABLE 3-1 
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TABLE 3-2 



LATERAL DIRECTIONAL STABILITY DERIVATIVES 
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CHAPTER 4 



LONGITUDINAL D Y NAM IC S 



4.1 INTRODUCTION 

Having specified all the vehicle's dynamic properties, we 
will now examine the dynamic motion, i.e. the vehicle's 
response to disturbances and control inputs. 

Since any flight vehicle is a dynamic system, described 
by differential equations of time, we must expect in a 
particular situation, oscillatory or exponential converging 
or diverging response. 

Any particular response will be associated by its 
corresponding characteristic root location in the complex 
plane, as we have clearly seen in the dynamic response of the 
spring mass damper system. 

The characteristic polynomials in the Laplace variable 
(s) involved, are usually of the third degree or higher and 
thus hard to solve by hand. 

In the past, a lot of approximations were introduced and 
even reduction of airframe's degrees of freedom technique was 
required to analyze coarsely a particular mode of motion, by 
reducing in effect the degree of the characteristic 
polynomial . 

In this chapter, we will analyze the longitudinal motion 
of a conventional airplane. We will first discuss the 
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Laplace Transform solution by obtaining the roots of the 
associated characteristic polynomial. Then we will apply the 
reduction of degrees of freedom technique to isolate and 
analyze a particular mode of motion. Finally, the equations 
of motion will be rearranged in a form that is more suitable 
for computer solutions. 

As an example, the longitudinal dynamic motion of the 
aircraft shown in Figure 4-1 will be analyzed with the 
parameters indicated. 

4.2 CHARACTERISTIC POLYNOMIAL OF LONGITUDINAL MOTION 

Longitudinal equations of motion are farther simplified 
by neglecting the effect of the Zq and Zw derivatives. 
Experimental results have shown that the effect of neglecting 
them is only of secondary importance. Further a straight 
level flight is assumed as the reference flight condition 
(i.e. =0) , hence. Equation 2-152 can be reduced to: 



S • Xy 0 Q 




V 




xa 


~ Zv s-Zw -ty> o 




w 




zs 


-Mu -sMw-Mw s-M^ 0 




% 




MS 


0 o i - s 




0 




0 



The characteristic polynomial is found by evaluation of 
the above determinant.' The Equation 4-1 can be reduced in 3 
equations of the original u,w,q, longitudinal variables by 
avoiding to write down the definition equation for q, but 
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attitude Cft) 20000 

awighl Clb) SOSOO 

macb number 0.638 

airspeed CWiec) €60 

X 0 = -0.0097 l, « -0.09SS Uu » 0.0 

X» = 0.0016 Zw* -< 45 *4* » -0.0235 

X$ * 0.0 Zb = -69.5 Mi* -0.0013 

-1.92 
MS 3 - 26.10 






( all dimensions in -fO 



Figure 4-1 

Conventional Airplane and Parameters 
Used for the Numerical Example 
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that doesn't simplify the hand calculations. The determinant 



comes out to be: 

As 4 + Bs 3 + Cs 2 + U s + £ =0 (4-2) 

where 

A-i (4-3) 

B = -Zvg -Xu (4-4) 

C = Zw - U,Mw ~ Xw Zv +• Xo ( M<^ + U,M w *• Zw) ( 4—5 ) 

U - Xo C - U> Mw) + Zu(Xw fej Mw) - Mv ( U>Xw ) ( 4—6) 

E = | (Z„M* - MuZ*) (4-7) 



Equation 4-2 is referred to as the longitudinal stability 
quartic. The roots in this case for any conventional vehicle 
come out to be complex with negative real parts, which are 
associated with two oscillatory convergent motions. 

The overall longitudinal motion of the vehicle is a 
superposition or a combination of these two oscillatory 
modes . 

One of the modes usually is associated with a pair of 
complex roots located near the origin of the negative complex 
plane. This mode is called "Phugoid motion" and is a fairly 
lightly damped oscillation with long periodic time. 

The other mode usually is associated with a pair of 
complex roots located far from the origin of the negative 
complex plane. This mode is called "Short period motion" and 
is a heavily damped oscillation with very short periodic 
time . 
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The longitudinal oscillatory modes are shown in 
Figure 4-2. 




The Longitudinal Oscillatory Modes 



For the given example we find 

A = 1.0 (4-8) 

B = 4.2177 ( 4_9 > 

C = 18.2962 (4-10) 

D = 0.1814 (4-11) 

E = 0.0722 (4-12) 

and the characteristic polynomial is 

S 4 + 4.2177 s 3 + *8. 29S2 s" + 0.1814 s f 0.0722 =0 (4-13) 

which has roots 

-0.0045 ± j 0.0627 (4-14) 

-2.1043 ± j 3.7184 (4-15) 
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The pair of roots -0.0045 ± j 0.0627 correspond to the 
phugoid motion and the pair of roots -2.1043 ± j 3.7184 
correspond to the short motion. 

The fourth degree polynomial can be factored in two 
second degree polynomials. Each one of the second degree 
polynomials is associated with the corresponding mode. 

Further it is more convenient to put each one of the 
polynomials in the familiar form 

5 2 + 2?u)ti S +■ W (4-16) 

so that the natural frequency c 0* , the damping 5 are figured 
out by inspection while the damped frequency can be 
determined by: 

(0 d = con ( (-5* y' Z (4-17) 

If we denote with the subscript (SP) the short period 
quantities and with the subscript (PH) the phugoid 
quantities, the characteristic polynomial may be written in 
the form: 

($ 2 + 25^ C0n sp S <-C0 n 2 sp )(s z t 25* S tCOn p4( ) =0 (4-18) 

Those quantities and their relation is shown in the 
complex plane in Figure 4-3. 

The roots of the characteristic polynomial are clearly 
shown in Figure 4-3 to be 

- u>rip H ± j tU df|4 (4-19) 

- 5 SP cur, sp i j (4-20) 
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or 

- 5 W U> n?H ± j 0)ri FH ( 1- 5* )' /2 ( 4-21 ) 

~ Ssp ^sp * j U)n S p C 4 ~ 5 ^ (4—22) 




Figure 4-3 

Longitudinal Roots and Characteristic Quantities 
in the Complex Plane 

Other quantities of importance in the longitudinal 
dynamics are: The phugoid period denoted by T vu , the short 

period period T SP given by: 





= 2rr/oj dM 


(sec) 


(4-23) 


Tsp 


= 2n/aj dsp 


Csec) 


(4-24) 



and the time for half amplitude for phugoid and short period 
given by: 
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0.6931 

5 p * (tin fM 



(4-25) 



0.6951 



-V 2 . 



:sP 5 spWt, S p (4-26) 

The Equations 4-25 and 4-26 come from the fact that the 
oscillatory response decays exponentially described by the 

-Ret 

envelope e where Re is the real part of the 

characteristic root. Therefore, to obtain half of the 
amplitude the relation 

<2 = 0.5 (4-27) 



must exist. By taking the natural logarithm of both sides we 
obtain the above relations. 

For our example, the characteristic polynomial of 
Equation 4-13 can be factored in the form 

(s 2 + 0.0090 s + 0 .00395) (s 2 + 4.2086s + 18.2546) (4-28) 



From this, by inspection we can determine: 

Wr> w = (0.00395) </2 = 0.0628 rad/sec (4-29) 

Wn sp = (18.2546) j/2 = 4.2725 rad/sec (4-30) 



$ 



pu 



0.0090 
2 * 0.0628 



0.0717 



(4-31) 




and also 



4.2086 

2x4.2725 



0.4925 



(4-32) 




0.6931 

0.0717 ' 0.OG28 



155-95 sec 



(4-33) 
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sec 



(4-34) 



= 0.0627 rad/sec (4-35) 

= 3.7184 rad/sec (4-36) 

(4-37) 

(4-38) 

It is worthwhile to note that: 

(a) The damped frequency is indicated by the 
imaginary part of the corresponding characteristic root. 

(b) The real part of the characteristic root 
indicates the exponentially decaying envelope of the 
corresponding oscillatory motion, i.e. the damping. 

(c) The following relations hold 



15 = 2 uJ,, + 2 5 sp u>n 5 p 




(4-39) 


C = (lJripy + A 5 S p cUr, sp 


+ 


(4-40) 


- 2 COn 4-2 


sp W r) p,j 


(4-41) 






(4-42) 



The two types of dynamical longitudinal oscillations are 
shown in Figure 4-4. 

The forward velocity v and the pitch angle © are usually 
excited in the phugoid mode while the w velocity and, hence, 
the angle of attack «• remains nearly unchanged. 






0,6951 

SP ~ 0.4925 * 4.2725 



= 0.55 



(% ?u = 0. 0628(1-0. 0717 2 ) 



2 ^Iz 



Cdfisp = 4. 2725(1-0. 4925 2 ) 



2 J/2 



2n 



PH = 



0.0627 



= -100.21 sec 



2n 



SP " 



3.71 84 



= -1.69 sec 
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Figure 4-4 

Dynamic Longitudinal Oscillations 

In contrast, the angle of attack is heavily excited in 
the short period mode as well as the pitch angle 0 , while the 
forward velocity remains nearly unchanged in that mode. 

The pilot can generally control the phugoid oscillation 
although because of his slow reaction it is possible to 
worsen to unstable due to out of phase action of his control. 
Phugoid oscillation must be highly considered in take-offs 
and landings. Short period is usually out of pilot control, 
but is of importance in automatic control. 

The requirements that an aircraft should satisfy are 
written in appropriate manuals. It is usually imposed that 
the oscillation should not take more than one cycle to damp 
to half amplitude. 



164 



Flight test can be carried out to determine the 
oscillations characteristics under different conditions of 
flight. Gyroscopic instruments can be placed and recordings 
of the longitudinal variables can be made possible to obtain. 
The aircraft is trimmed to the selected speed and height and 
longitudinal oscillation is initiated by an abrupt elevator 
deflection. 

4.3 PHUGOID AND SHORT PERIOD APPROXIMATIONS 

The fact that the angle of attack and hence, the 
w-velocity remains nearly unchanged in the phugoid mode, gave 
the idea of determining the phugoid quadratic by: 



(a) eliminating the w terms. 

(b) eliminating the pitching moment equation since no 



considerable pitching moment variation occurs, and 

(c) letting q=s9, obtaining thus the determinant: 



-Zo -UoS 



(4-43) 



The characteristic phugoid quadratic comes out to be: 




'O 



(4-44) 



For example, given, the characteristic quadratic is 



s 2 +0. 0097s + 0 .0047 = 0 



(4-45) 



which has roots 



-0.0049 ± j 0.0681 
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(4-46) 



These roots are approximately equal to the exact phugoid 
roots . 

From Equation 4-44 it is clear that the natural frequency 
of the phugoid is heavily dependent upon the Zu derivative 
while the damping of the phugoid is dependent upon the X u 
derivative. It is also obvious that high speeds tend to 
reduce the phugoid frequency and to increase the 
corresponding period. 

Since for subsonic straight level flight the Zu 
derivative is simply expressed as: 

Zo=-2%lDo (4-47) 

the following approximate relations hold for the phugoid 







(4-48) 



t Do Xu 

bpW = ~ Q\l? £ (4-49) 

Aircraft designers, in order to achieve larger values of 
phugoid damping, must increase the value of the Xu derivative 
but this derivative is proportional to the drag. 
Fortunately, in take-offs and landings where phugoid motion 
is undesired, high drag devices such as flaps are used 
anyway . 

A more precise approximation for the phugoid quadratic is 
given by leaving the w terms and the Hu term which in modern 
crafts is never zero. This is a three degree of freedom 
phugoid approximation. 
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s-Vo y w ^ 

-Zo s-Zvu _U 0 (4-50) 

- Mo -Mw 0 

The characteristic phugoid quadratic becomes then 
Mo (UoXw -j) 1 g f Mo 



s 2 + 



f y (.UoXw -|) *1 g ( Mo - \ 



(4-51) 



which gives 



= ri- 1 ^ z .,)]" 1 



and 



5 j 



[~Xt> 



M*/ 

Mo (L) 0 Xw -a) 
UoMw 



(4-52) 



■j 



>[£< 4 -anr 



(4-53) 



Short period approximations can be made possible by: 

(a) eliminating the v-equations since very small 
v-velocity change occurs. 

(b) setting all X-force derivatives to zero since 
very small X-force variation occurs. 

(c) letting q=s9, obtaining thus the determinant: 

■ s~ Zw ~ Do 

-sMw-Mw s-M^ (4-54) 

The characteristic short period quadratic comes out to 



be : 



s 2 + (-U>Mw - -Z w ) s + (-U,Mw + M^Zvu) = 0 (4-55 ) 
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